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^0.  Abstract  (continued) 

treated  independently.  Each  impressed  current  produces  a  "dominant"  current 
distribution,  a  characteristic  of  the  element,  which  can  be  readily  computed. 
Current  coupling  is  formulated  by  "intrinsic"  impedance  matrices  which  relate 
the  scalar  potentials  at  the  terminals  of  an  element,  caused  by  its  dominant 
current  distributions,  to  the  impressed  currents  of  the  element.  Field  coupling 
produces  "scatter"  currents  on  all  the  elements,  and  is  formulated  by  a  "field 
coupling"  matrix  which  relates  the  scalar  potentials  at  the  terminals,  caused 
by  field  coupling,  to  the  impressed  currents  at  all  the  terminals.  Intrinsic 
and  "field  coupling"  are  combined  to  form  the  "complete"  impedance  matrix  of 
the  diakopted  antenna.  Enforcing  continuity  of  the  currents  and  equality  of 
the  scalar  potentials  at  all  the  interconnections  between  the  elements  yields  a 
system  of  linear  equations  for  the  junction  currents  and  the  input  impedance  of 
the  antenna.  Current  coupling  dominates  over  field  coupling.  Field  coupling 
due  to  the  dominant  current  distributions  of  the  elements  is  of  primary  importance 
wnile  field  coupling  due  to  the  scatter  currents  is,  in  general,  negligible. 

This  theory  is  applied  to  several  multi-element  antennas  and  the  results  are 
compared  with  other  methods  to  highlight  the  numerical  advantages. 
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Summary 


This  research  document  presents  a  new  theory  for  the  analysis  of  multi¬ 
element  antennas  which  consist  of  interconnected  conductive  structure  elements 
of  electrically  small  dimensions.  The  theory  is  based  on  the  retarded  electro¬ 
magnetic  potentials  which  permit  a  diakoptic  approach  to  the  problem.  The 
antenna  is  broken  up  into  its  individual  structure  elements.  Each  element  is 
assumed  to  be  excited,  a)  by  currents  which  are  impressed  at  its  terminals, 
i.e.  junctions  with  adjacent  elements  (current  coupling),  and  b)  by  the  electric 
fields  of  the  currents  and  charges  on  all  the  other  elements  (field  coupling). 
Both  excitations  are  treated  independently.  Each  impressed  current  produces  a 
"dominant"  current  distribution,  a  characteri sti c  of  the  element,  which  can  be 
readily  computed.  Current  coupling  is  formulated  by  "intrinsic"  impedance 
matrices  which  relate  the  scalar  potentials  at  the  terminals  of  an  element, 
caused  by  its  dominant  current  distributions,  to  the  impressed  currents  of  the 
element.  Field  coupling  produces  "scatter"  currents  on  all  the  elements,  and 
is  formulated  by  a  “field  coupling"  matrix  which  relates  the  scalar  potentials 
at  the  terminals,  caused  by  field  coupling,  to  the  impressed  currents  at  all 
the  terminals.  Intrinsic  and  "field  coupling"  are  combined  to  form  the 
"complete"  impedance  matrix  of  the  diakopted  antenna.  Enforcing  continuity  of 
the  currents  and  equality  of  the  scalar  potentials  at  all  the  interconnections 
between  the  elements  yields  a  system  of  linear  equations  for  the  junction 
currents  and  the  input  impedance  of  the  antenna.  Current  coupling  dominates 
over  field  coupling.  Field  coupling  due  to  the  dominant  current  distributions 
of  the  elements  is  of  primary  importance  while  field  coupling  due  to  the 
scatter  currents  is,  in  general,  negligible.  This  theory  is  applied  to  several 
mul ti -element  antennas  and  the  results  are  compared  with  otner  methods  to 
highlight  the  numerical  advantages. 
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I.  Introduction 


Improved  tactical  communication  systems  require  antennas  which  are 
electrically  small  (i.e.  small  compared  witf,  the  wavelength),  have  very  large 
bandwidths  and  reasonably  high  efficiency.  It  is  well  known  to  antenna 
experts  that  these  requirements  work  against  each  other.  The  problem  therefore, 
is  to  find  sophisticated  antenna  structures  which  provide  the  best  compromise 
between  these  contradicting  requirements. 

Experimental  investigations  of  empirically  designed  multielement  antennas, 
i.e.,  antennas  which  comprise  a  number  of  interconnected  and  closely  spaced 
conductive  elements,  have  shown  promising  results.  An  example  of  such  a  broad¬ 
band  multielement  monopole  antenna  is  shown  in  Figure  1.  This  antenna  con¬ 
sists  of  four  vertical  conductors.  The  two  thicker  ones  are  grounded,  while 
the  other  two  are  interconnected  near  the  ground  plane  and  connected  to  the 
input  terminal.  Each  vertical  conductor  has  a  top  capacitor  in  the  form  of 
a  metal  plate,  and  there  are  inductive  interconnections  between  the  plates  in 
the  form  of  wire  loops.  But  antennas  like  the  one  mentioned,  whose  functioning 
is  not  quite  understood,  are  not  amenable  to  conventional  computer  analysis. 

An  analytical  treatment  of  such  a  composite  structure  appears  to  be  a 
rather  hopeless  undertaking.  Commonly  used  numerical  techniques  are  impracti¬ 
cable  because  they  would  require  computers  with  enormous  storage.  Moreover, 
these  techniques  do  not  always  yield  reliable  results  [2]. 

This  research  offers  a  new  approach  to  problems  of  this  kind.  According 
to  this  approach  the  composite  structure  is  diakopted  into  its  individual 
structure  elements.  As  a  simple  example,  Figure  2  shows  a  diakopted  dipole 
with  end  capacitor  plates.  Each  structure  element  is  characterized  by 
electrical  quantities  which  depend  only  on  size  and  shape  of  the  element,  and 
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the  assembly  is  treated  similarly  to  the  interconnection  of  n-port  networks. 

The  excitation  of  each  element  is  ascribed  to  two  causes,  a)  the  currents 
entering  the  element  at  its  "terminals,"  i.e.  junctions  with  adjacent  elements 
or  the  source,  and  b)  the  fields  of  the  currents  and  charges  on  all  the  other 
elements.  The  first  is  referred  to  as  "current  coupling"  and  the  second  as 
"field  coupling."  Both  excitations  are  treated  separately.  Current  coupling 
implies  hypothetical  sources  with  a  single  terminal  and  the  capability  of 
impressing  a  current  onto  a  conductor.  Although  such  sources  violate  the  conti¬ 
nuity  condition,  their  assumption  is  permissible  if  the  electro-magnetic  fields 
are  expressed  by  the  retarded  electromagnetic  potentials.  Although  the  conti¬ 
nuity  condition  is  violated  in  the  treatment  of  individual  structure  elements, 
it  is  restored  when  the  elements  are  interconnected.  Thus,  current  coupling  is 
computed  by  impressing  a  current  at  a  terminal  of  a  structure  element.  This 

current  spreads  over  the  surface  of  the  element  and  produces  a  current  distri¬ 
bution  which  is  uniquely  determined  by  the  geometry  of  the  element  and  the 

location  of  the  terminal  and  is  called  the  dominant  current  distribution 

associated  with  a  given  terminal.  There  are  as  many  dominant  current  distri¬ 
butions  as  there  are  terminals.  The  relationship  between  the  scalar  potentials 
at  the  terminals  (produced  by  the  dominant  current  distributions)  and  the  im¬ 
pressed  currents  is  formulated  by  the  "intrinsic  impedance  matrix"  of  the 
element. 

Field  coupling, on  the  other  hand, excites  scatter  currents  which  are  super¬ 
imposed  on  the  dominant  current  distributions.  The  scalar  potentials  at  the 
terminals  due  to  field  coupling  depend  on  all  the  impressed  currents.  Their 
relationship  with  the  impressed  currents  is  forumlated  by  a  "field  coupling" 
matrix.  The  intrinsic  impedance  matrix  and  "field  coupling"  matrix  combined 
together  form  the  "complete  impedance  matrix"  of  the  diakopted  antenna.  This 
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matrix  relates  the  total  scalar  potentials  at  the  terminals  of  all  the  elements 


to  all  the  impressed  currents. 

Interconnection  of  the  structure  elements,  which  requires  equal  scalar 
potentials  at  the  interconnected  terminals  and  continuity  of  the  junction 
currents,  it  formulated  by  an  interconnection  matrix.  In  this  manner  a  system 
of  linear  equations  is  obtained  which  yields  the  junction  currents  and  the 
input  impedance  of  the  antenna. 

A  most  simple  antenna  to  which  the  theory  aDDlies  is  a  simple  monopole 
antenna  with  a  top  capacitor.  In  this  case,  there  are  two  structure  elements, 
the  vertical  conductor  and  the  top  capacitor.  The  qround  plane  can  be  replaced 
by  the  antenna  image.  No  systematic  way  of  computing  the  impedance  character¬ 
istics  of  this  antenna  has  been  reported  in  the  literature. 
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1 1 .  Diakoptic  Theory  of  Mul ti -Element  Antennas 

In  this  section  we  shall  develop  the  essential  theoretical  results  re¬ 
quired  to  implement  the  diakoptic  theory. 

Consider  a  mul ti-element  radiating  structure  such  as  shown  in  Figure  1. 
Various  elements  are  interconnected  to  each  other  via  terminals  of  junctions. 
Let  each  radiating  element  be  disconnected  or  (diakopted)  from  all  other  ele¬ 
ments  and  be  suspended  in  space.  The  assemblage  of  these  disconnected  ele¬ 
ments  is  called  the  diakopted  (or  primitive)  system.  Each  element  has  many 
terminals  on  each  of  which  certain  impressed  current  and  potential  is  assumed. 
The  essential  requirement  for  this  diakopted  system  with  impressed  currents 
along  the  junctions  is  that  it  be  performancewi se  identical  to  the  assembled 
antenna.  Thus, 

a)  The  sum  of  the  impressed  currents  is  zero  at  every  junction  between 
the  structure  elements  and  the  continuity  condition  is  satisfied  at 
every  input  terminal.  This  requirement  assures  that  the  field  of 
assembled  antenna  is  Maxwellian. 

b)  The  scalar  potentials  at  the  interconnected  terminals  are  equal. 

c)  The  potential  difference  between  the  input  in  terminals  is  equated 
with  the  driving  voltage  of  the  antenna  source. 

Let  the  potential -current  relationship  at  every  terminal  be  written 


in  matrix  form: 

[q]  =  [Z ] C 1 3  [diakopted  antenna]  II. 1 

[$]’  =  [Z 3 ' [ I ] '  [assembled  actual  antenna]  II. 2 

Requirements  (a),  (b)  and  (c)  represent  Kirchoff's  laws  for  interconnected 
structures  and  can  be  written  as 

[1]  =  [C][I]‘  •  --  II. 3 

c£r  =  cc]t[$]  ii. ^ 
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II. 5 


[,']'tci]'  -  ci]tm 

[C]  represents  the  transpose  [C], 

Matrix  [Z]  represents  the  impedance  of  the  diakopted  antenna  and  primed  quantities 
refer  to  the  actual  assembled  antenna.  [C]  may  be  a  rectangular  matrix  with 
(C.j  .}  os  0  or  1. 

From  II. 3,  II. 4  and  II. 5,  the  impedance  matrix  of  the  actual  structure  can  be 
written  as 

[Z]1  =  [CJt[Z][C] 

An  example  at  the  end  of  this  section  shows  how  [C]  and  [Z J '  are  obtained. 

The  essential  results  of  this  section  show  that  in  order  to  obtain  [Z]' , 
we  have  to  only  compute  the  impedance  matrix  [Z]  of  the  so  called  diakopted 
structure. 

The  most  important  point  here  to  remember  is  that  the  elements  of  the 
impedance  matrix  [Z]'  depend  upon  simultaneously  knowing  current  distribution 
on  all  the  radiating  structure  elements.  Thus,  without  diakopting  the  struc¬ 
ture,  we  have  to  simultaneously  solve  as  many  integral  equals  as  there  are 
radiating  elements.  On  the  other  hand,  the  elements  of  the  impedance  matrix 
[Z]  of  the  diakopted  structure  can  be  found  by  computing  the  current  distribu¬ 
tion  on  individual  elements  separately  and  hence  involves  solving  as  many 
integral  equations  as  there  are  radiating  structures,  but  only  individually. 

This  results  in  a  tremendous  savings  of  numerical  computation.  In  what  follows 
we  shall  show  how  the  so  called  total,  primitive  (or  diakopted)  impedance 
matrix  [Z]  can  be  computed. 
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III.  Impedance  Matrix  of  a  Diakopted  Antenna 

Consider  each  element  of  a  diakopted  antenna. 

The  excitation  of  each  element  is  ascribed  to  two  causes,  a)  the  currents 
entering  the  element  at  its  "terminals,"  i.e.,  junctions  with  adjacent  elements 
or  the  source,  and  b)  the  fields  of  the  currents  and  charges  on  all  the  other 
elements.  The  first  is  referred  to  as  "current  coupling"  and  the  second  as 
"field  coupling."  Both  excitations  can  be  treated  separately  and  the  resulting 
compling  can  be  superimposed  due  to  linearity.  Current  coupling  implies  hypo¬ 
thetical  sources  with  a  single  terminal  and  the  capability  of  impressing  a 
current  onto  a  conductor.  Although  such  sources  violate  the  continuity  condition, 
their  assumption  is  permissible  if  the  electro-magnetic  fields  are  expressed  by 
the  retarded  electromagnetic  potentials.  Although  the  continuity  condition  is 
violated  in  the  treatment  of  individual  structure  elements,  it  is  restored  when 
the  elements  are  interconnected.  If  a  current  is  impressed  at  a  terminal  of  a 
structure  element, the  current  spreads  over  the  surface  of  the  element  and  pro¬ 
duces  a  current  distribution  which  is  uniquely  determined  by  the  geometry  of  the 
element  and  the  location  of  the  terminal.  This  is  called  dominant  current  dis¬ 
tribution  of  a  particular  element.  Thereareas  many  dominant  current  distribu¬ 
tions  as  there  are  terminals.  The  relationship  between  the  scalar  potentials  at 
the  terminals,  produced  by  the  dominant  current  distributions,  and  the  impressed 
currents  is  formulated  by  the  "intrinsic  impedance  matrix"  of  the  element  and 
is  referred  to  as  [Z(I)j. 

Field  coupling  excites  scatter  currents  which  are  superimposed  on  the  domi¬ 
nant  current  distributions.  The  scalar  potentials  at  the  terminals  due  to  field 
coupling  depend  on  all  the  impressed  currents.  Their  relationship  with  the  im¬ 
pressed  currents  is  formulated  by  a  "field  coupling"  matrix  [Z(F)]. 

[Zj  *  [Z(I)]  +  [Z(F)j 


10 
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This  matrix  [2]  is  called  the  total  impedance  matrix  of  the  diakopted  antenna 
and  relates  the  total  scalar  potentials  at  all  the  terminals  of  all  the  ele¬ 
ments  of  the  diakopted  structure  to  all  the  impressed  currents. 

1 1 1  - 1  Current  Coupling  Between  Structure  Elements  and  Intrinsic 
Impedance  Matrix  [Z(I)]. 

A.  Structure  elements  with  one  terminal 

Consider  one  of  the  capacitor  plates  of  the  dipole  in  Fig.  2  separated 
from  the  other  elements  and  suspended  in  space,  with  a  current  I  impressed  at 
the  terminal,  i.e.,  contact  area  in  the  center  of  the  plate  (Fig.  3).  The  con¬ 
tact  area  a  is  considered  very  small  compared  with  the  surface  area  of  the 
element.  Excitation  by  an  impressed  current  cannot  be  treated  with  Maxwell's 
equations,  because  Maxwell's  equations  imply  sources  which  separate  positive 
and  negative  charges.  In  contrast,  impressed  currents  require  sources  which 
produce  charges.  The  retarded  electromagnetic  potentials  do  not  impose  any 
conditions  on  the  source,  and  can  therefore  be  used  for  our  problem. 

If  i(r)  is  the  surface  current  density,  and  q(r)  the  surface  charge 
density  due  to  the  impressed  current  I,  the  retarded  potentials  are 

A(r)  =  1  T(? ' )G(? ,? ' )dS  (vector  potential)  III.l 


:<r) 


-. —  q(r')G(r,r')dS  (scalar  ootential) 

■Xl,e  J  S 


III. 2 


wi  th 


j(r,r ' ) 


exp(  - jk 1 r  -  r ' i ) 


k  =  2-7  x 


! 


h 

i  \ 


wntre  r'  is  the  position  vector'  or  the  chorees  uric  cur-'ents  fie 
elements  dS ,  and  r  that  of  the  point  of  observation.  The  quantities  i  ro 
q(r)  must  satisfy  the  following  two  eauations  on  the  surface  of  toe  e ent 
outside  the  contact  area  : 

E(r)  ' d?  =-  [jwA(f)  +  7a(r)]'  dS  =  0  (Boundary  condition,'  III. 3 

7-  T( r )  +  juiq(r)  =  0  (Continuity  condition'  III. 4 

The  condition  that  the  current  flux  through  the  boundary  curve  7  of  t.ne  con¬ 
tact  area  o  is  the  continuation  of  the  impressed  current  I,  is  given  as: 

C  T(r) *t(r)dr  =  I  III. 5 

where  t(r)  is  a  unit  vector  tangential  to  the  surface  S  and  normal  :.  ~he 
current  and  charge  distributions  T(r)  and  q(r)  due  to  the  impressed  current  I 
are  termed  as  "dominant"  distributions  since  the  currents  due  to  field  coupling 
between  the  elements  are,  in  general,  relatively  small.  From  the  boundary 

condition  III. 3 

!  E(r)*?(r)dS  =  -  I  [jwA(r)  +  7$(r)]*?(r)dS  =  0  III. 6 

'  S  J  S 

The  surface  of  integration  S  is  the  surface  of  the  element  with  the  exclusion 
of  the  contact  area.  Using  the  relations 

7p  •  i  ( r)  =  7-(o(r)!(r) )  -  $(r)v-?(r)  =  V*(>(r)!(r);  +  ><i(r)q(r)  IH-7 

and  applying  Gauss'  theorem,  one  obtains  from  III. 6 


i.A(r)-?(r}  +  ?(r)q(r)]d  $  =  -  1  7*  ($(  r}?{  r)  }cs  =  O  :  (  r)  i ;  r)  •  t(  r]  dr  III. 8 
•  S  •  S 

I*  the  contact  area  .7  is  sufficiently  small,  0  can  be  considered  constant  w  itn- 
;n  the  contact  area.  Thus,  with  (5),  equation  (3)  reduces  to 


juf  [A(r)*T(r)  +  o(r)q(r)ldS  =  $1 
'S 


III. 9 


w.nere  >  is  the  scalar  ootential  at  the  contact  area. 

The  ratio  between  :>  and  I  can  be  used  to  define  an  impedance  which  shall 
oe  termed  "intrinsic  impedance."  If  A  and  t  are  expressed  by  tne  current  and 
z-'arge  distribution,  the  intrinsic  impedance  of  the  element  is 


Z(r)  =  f  =  ij  (  [A( r) »T( r)  +  i>(r)q(r)]dS 
1  r  h 


«^f  f  G(?,r 1 )  _  -L  SlflolLll 


JS  JS‘ 


Qc 


dS 


III. 10 


where  Q  =  I/joj  is  the  total  charge  on  the  element.  The  current  and  charge  dis¬ 
tribution  functions  i / I  and  a/Q  are  solely  determined  by  the  geometry  of  tne 
element  and  the  location  of  the  coupling  area. 

When  the  intrinsic  impedance  is  computed  with  Ili.io  for  a  conductor  of  any 
shape,  for  extremely  low  frequencies,  it  takes  the  form 

Z(I)  •  jCc  •  47/F  <“-0>  III.ll 


where  C  is  the  static  capacity  of  the  element.  The  first  term  1/j^C  is  the  one 
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wmcn  is  to  be  expected.  The  second  tern  represents  a  negative  res: stance 
-lj  onms  and  is  not  quite  obvious.  It  is  brought  about  by  the  ’"‘act  t.-.et  an 
impressed  current  produces  a  charge  on  tne  element  without  a  counterc-.arge . 
in  contrast  to  a  Maxwell  source.  If  the  scalar  potential  is  expanded  in  a 
oower  series  in  u,  one  obtains 


.(r 


j  G(r,r')q(r‘)dS'  = 
JS' 


V 


jk'  q(r'}cS' 

;  C  1 
o 


The  first  term  of  this  expansion  is  the  static  potential  of  the  cnarges.  The 
second  term  which  is  independent  of  r  represents  a  potential,  termed  "bacx- 
ground"  potential  $  ,  which  is  uniform  in  space  and  has  no  gradient.  ~his 
means  it  does  not  produce  a  field.  It  is  this  background  potential  which 
produces  the  -30ft  term  in  III.  11.  When  the  element  which  we  assumed  to  be  sus¬ 
pended  in  space  is  within  the  antenna  structure  the  background  potential  is 
compensated  because  the  combined  charges  on  all  the  other  elements  are  negative¬ 
ly  equal  to  the  charges  of  the  considered  element.  The  background  potential 
can  be  avoided  if  the  retarded  scalar  potential  is  redefined  as  modified  scalar 
potential 


where 


J_  f 

4ne  Js 


G(r,r ' )q(r 1 )dS  ' 


G(r,r' ) 


-jk|r-r' ! 

—  .  ■  +  Jk 

!  v*  -  ^  1  ! 


HI.  ;2 


This  modified  scalar  potential  which  will  be  used  throughout  the  pacer  -is  legit¬ 
imate  as  it  is  not  conflicting  with  Maxwell's  theory.  Since  =  ?:•,  the 


boundary  condition  III.  3  ana  the  dominant  current  distribution  cer-  -  : 
tnerer'rotn  remain  unchanged  if  tne  conventional  potential  :  is  substi t-teo 
the  modified  potential  ;..For  a  Maxwell  system  i  and  ;  are  identical,  since  ccC 
extended  over  the  entire  surface  of  the  system  is  zero.  The  intrinsic  ir.oedance 
of  a  structure  element  with  one  connection  becomes 


7  = 


[A (r)-i(r) 


+  f(r)q(r)]dS 


]UU 

4tt 


f  I  !G(?,r‘ 

jSJS  ! 


)  Umtl  .  1  G(r,r 1 ) 

r  k~  o4 


dG  ‘  dS 
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III. 14  represents  a  stationary  formulation  of  the  intrinsic  imped¬ 
ance.  This  means,  small  errors  in  the  dominant  current  distribution  have  only 
a  second  order  effect  on  the  intrinsic  impedance.  (See  Appendix  J) 

Excitation  by  an  impressed  current  I  at  the  terminal  can  be  considered 
equivalent  with  the  excitation  by  an  oscillating  charge 


Q 


jw 


III. 15 


which  is  placed  above  the  contact  area  at  a  distance  d^C  as  shown  in  Fig.  a. 

The  charge  on  the  contact  area  a  consists  essentially  of  the  image  charge  -Q, 
with  the  charge  +Q-distributed  over  the  surface  areas  of  the  structure  element, 
because  the  net  charge  on  the  element  must  be  zero.  The  equivalence  between 
charge  and  current  excitation  is  shown  in  Appendix  1. 

The  intrinsic  impedance  Z(I)  of  an  element  with  one  terminal  can  be  repre¬ 
sented  by  a  lumped  element  circuit  as  snown  in  Fig.  5.  For  low  frequencies, 
•.e.  wnen  the  dimensions  of  the  element  are  small  compared  with  tne  wave  ' 

:  and  L  can  be  considered  constant,  while  R  increases  proportionally  with  oj2: 

Z  =  tV  +  +  R(w“} 


,*7»-  . 


In  case  *p  instead  of  p  is  used,  the  individual  elements  will  show  an  additional  -3011. 
fictitious  resistance  in  the  intrinsic  impedance.  However,  the  impedance  of  the 
totally  assembled  antenna  is  the  same  as  the  conventional  impedance  due  to 
automatic  compensation  of  -30ft. 


3 .  Structure  elements  with  two  or  more  terminals 

A  structure  element  with  two  terminals  such  as  th^  cylindrical  conductors 
of  Fig.  2  has  two  dominant  current  distributions,  one  associated  with  each  of 
the  impressed  currents  (Fig.  6).  Each  dominant  current  distribution  oroduces 
a  scalar  potential  at  both  contact  areas.  If  t.-j  and  are  the  potentials 
at  the  terminals  1  and  2  due  to  1  ^ ,  and  those  due  to  I9,  then,  the 

relationship  between  the  total  potentials  4-j  and  4^  at  the  terminals  and  the 
impressed  currents  can  be  written  as 


*1  =  *11  +  *12  =  Z11(I)I1  +  Z12^I)I2 
*2  =  *21  +  *22  =  Z12^I^I1  +  Z22^*2 


iii.i; 


For  a  structure  element  with  M  terminals  the  relationship  between  the 
terminal  potentials  and  the  impressed  currents  is  formulated  by  an  M  x  M  intrin¬ 
sic  impedance  matrix. 

CO  =  CZJCI]  III.  18 


wr.ere 


^k(r)qj-(r)]dS 
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J  SJ  S'  '  Yk  k2 


f  f 


q  ^  ( h )  q  k  (  r 1 )" 


dSdS 


III. 19 


18 


The  quantities  i.,q.  and  i,  ,q,  are  the  dominant  current  and  charge  distribu- 
J  J  K  K 

tions  generated  by  the  impressed  currents  I ^  =  jtoQj  and  1^  =  and  A^,^ 

are  the  retarded  potentials  associated  with  ik>q^.  (III. 19)  is  derived  in 
Appendix  2. 

The  symmetry  of  the  intrinsic  impedance  matrix,  Zjk(I)  =  Zkj(I),  is  evident 
from  the  second  formulation  of  III. 19.  In  Appendix  4  it  is  shown  that  III. 19 
is  a  stationary  representation  of  the  matrix  elements. 

A  lumped  element  equivalent  circuit  for  a  structure  element  with  two  ter¬ 
minals  is  shown  in  Fig.  7.  For  sufficiently  low  frequencies  the  capacitors  and 
inductors  can  be  considered  constant,  while  the  resistors  increase  with  oj2.  The 
resistor  which  is  in  series  with  the  capacitor  is  negative,  but  smaller  than  the 
resistors  associated  with  the  inductors. 

1 1 1 . 2  Field  Coupling  Between  Structural  Elements. 

Field  Coupling  Impedance  Matrix 

We  now  consider  a  diakopted  structure  and  arbitrary  currents  impressed  at 
the  terminals.  The  capacitively  loaded  dipole  of  Figure  2  may  serve  as  an  example. 
The  terminals  are  identified  by  a  superscript  i  and  a  subscript  k,  the  superscript 
referring  to  the  number  of  the  element,  and  the  subscript  referring  to  the  number 
of  the  terminal  on  the  element.  If  there  were  no  field  coupling  between  the 
elements,  the  current  distributions  on  all  the  elements  would  be  the  dominant 
distributions  associated  with  the  impressed  currents. 


The  field  of  a  dominant  current  distribution  is  non-Maxwell ian  since  the 
associated  net  charge  is  nonzero.  If  a  current  1^  is  impressed  at  the  terminal 
(k),  the  non-Maxwell ian  field  of  the  dominant  current  and  charge  distribution 
i^,  q^  induces  currents  on  all  the  other  elements.  The  scatter  fields  excited 
by  these  induced  currents  are  Maxwellian,  since  induced  current  distributions 
have  no  net  charge.  These  "first  order"  scatter  fields  excite  second  order 
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scatter  fields  and  so  on,  each  higher  order  having  a  greatly  reduced  amplitude. 

All  these  scatter  fields,  when  summed  up,  form  a  multiple  scatter  field  which 
is  Maxwellian.  The  currents  and  charges  associated  with  the  multiple  scatter 
field  are  distributed  over  all  the  surfaces  Sn  (including  S1 )  and  shall  be  denoted 
<5i?n,  q^.n,  the  super  and  subscripts  indicating  that  they  are  produced  by  the 
impressed  current  1^  and  located  on  tne  element  n. 

The  total  field  generated  by  i’  satisfies  on  every  element  the  bour.aar. 
condi tions 

(jw(A^  +  <sA^)  +  +  a<ji^))  x  dSn  =  0  (n  *  1,  ....  i,  ...N)  III. 20 

where  A^,  are  the  retarded  potentials  of  the  dominant  current  and  charge 
distribution  q’,  and  oA’,  5$’  those  of  the  scatter  current  and  charge  ais- 
tributions  ,  <sq^n  combined.  N  is  the  number  of  elements. 

Since  the  electric  field  of  T^,  q^  satisfies  the  boundary  condition  on 
S1,  it  follows  from  (20)  for  n  =  i  that 

(jw<5A^  +  V6$’)  x  dS1  =  0 

Thus 

!  ,•  |  ( jwoAji  +  75$?.)  •??  j  dS1  =  0,  '  III. 21 

^sM  k  k  kt  k  ■  1,  ....  M. 

' —  I 

Using  the  relations  III. 7  and  Gauss'  theorem  one  obtains  the  "backscatter" 
potential  due  to  the  field  interaction  of  the  excited  element  with  the 
other  elements: 
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>,(n  lk  ju,Li  (':Ak‘Tk+  5ikqk)dsl 

K  )  K  o 


in. 2; 


The  letter  F  indicates  field  coupling,  the  first  pair  of  indices  refers 
to  the  terminal  at  which  $  is  determined,  and  the  second  pair  to  the  terminal 
of  the  impressed  current  which  produces  this  potential. 

As  shown  in  Appendix  3 


.iuKr\ +  qk)ds1  =  X  U  (Ak*dk‘ +  vv 


.^tn  ,  • i .  in,cn 
(A.  *<n  *  +  •>.  -„q.  ;dS 


III. 23 


n=l  -'S' 


furthermore,  from  the  boundary  conditions  III. 20  using  the  relations  III. 7  and 

Gauss'  theorem,  follows 
r 

1  (A?  +  5 +  (■?/  +  5*?)  cqV’ldS0  =  0,  for  every  n  including  i 

j  ^n  •  K  X  K  K  K  K  | 


111.24 


The  right-hand  side  of  III. 24  is  zero  since  for  scatter  currents  tee 


integral  of  Gauss'  theorem  in  1 1 1. 5  is  zero  ( see  Appendix  5). 

From  the  last  three  equations,  one  obtains  the  "back  scatter"  impedances 


„  i  ,i 
i  i  «(F) 
Z(F)  =  -M-  = 


- 

l  k  ) 


1 1 1 . 25 


which  has  to  be  added  to  the  diagonal  terms  of  the  intrinsic  impedance  matrix 
Z ?  * ? •  using  the  notation  of  this  section.  Generalization  of  III. 25  to  obtain 

K  *  K 

the  scatter  field  contributions  to  the  off-diagonal  terms  is  straightforward. 
One  obtains 
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»-j*-TT  I  j  n(5Rk-di,n  +  3$j6q!’")dSn  Hi. 26 

1,1.  n=l  JSn  K  J  k  J 


For  k  =  j  equation  (26)  transforms  into  III. 25 


Let  us  now  determine  the  potential  $(F)  produced  by  the  impressed  cur- 

k,m 

rent  I2,  at  the  terminal  (l1).  Because  of  the  boundary  condition  1 1 1. 20 
m  k 


1  .[jw(A  +  6A  )  +  7(?  +  5<r)].T‘dS  =  0, 

j  2 1  m  m'  wm  m  k 


III. 27 


Cj«(Aj  +  oAj)  +  ml  +  5^)].6i;*ndSn  =  0 


III. 28 


where  the  second  equation  holds  for  every  n  including  i.  The  ootentials 


iA^  and  5^  characterize  the  scatter  field  which  would  be  excited  by  1^. 

As  before,  we  apply  III.  7  and  Gauss'  theorem  to  the  above  two  equations 


to  obtain 


5  '(•*-  La  +  <v>>sf  +  L(aK  *  sv(>  as,j 

k  ,m  i ;  5  5 


III. 29 


0  -  [sn<Si-5»n  *  »M>Sn  * 
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The  first  term  in  III. 29  represents  the  contribution  to  the  terminal  potential 
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i  ,<• 

;{F)  from  the  non-Maxwel 1 i an  field  of  the  dominant  current  and  cnarce  distribu- 
k  ,m 

tion  T‘  ,  q\  and  the  second  term  that  from  the  scatter  current  ana  c.narge  ais- 
m  m 

,  . .  ..  .ti,n  ~  f. , n 

tributions  ,  oq  ’ 

m  m 

As  shown  in  Appendix  3 
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J  ,A 

Expressing  $(F)  in  terms  of  an  impedance 
k,m 


i  ,1 


♦(F)  =  Z(F)  I 
k,m  k,m 


III. 32 


i  £ 

the  field  coupling  impedance  between  the  terminals  (^)  and  (  ) 
becomes 


i  ,l 
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r, 


Z(F)  =  i  * 
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1 1 1. 33 


Equations  III. 26  and  III. 27  formulate  the  elements  of  the  field  coupling  im- 

i 

pedance  matrix  [Z(F)]  which  relates  the  scalar  potentials  $(F)  at  the  terminals, 

k 

caused  by  field  coupling,  to  the  impressed  currents: 


M  M 

1  N  l  N  \  i,i  , 

[5(F)]  *  [Z(F)][I],  Hf)  =  l  l  9(f)  .  I  l  Z(F)  I  III. 34 

k  i* 1  m=l  k,m  i=i  m=i  k,m  m 
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The  intrinsic  impedance  matrices  of  the  individual  structure  elements 

can  be  combined  into  diagonal  block  impedance  matrix  [ Z ( 1 ) ]  by  writing  the 

matrix  elements  Z.  .  (Eq.  III. 19)  in  the  form  The  superscript  i 

K » J  k ,  j 

identifies  the  terminals  k  and  j  as  belonging  to  the  element.  The  block 

i  £ 

matrix  [Z ( I ) ]  whose  elements  Z. ’.(I)  are  zero  for  i  t  l  is  the  "current 

K  j  J 

coupling  matrix"  of  the  diakopted  system,  and  relates  the  terminal  potential 
N  M' 

♦Id)  =  i  i  4’i(I)  i\ 

K  i=l  j=l  K,J  J 

due  to  current  coupling  to  the  M.  impressed  currents  of  the  element  i. 

The  sum  of  the  matrices  [Z(I)]  and  [Z(F)],  i.e. 

[Z]  =  CZ(I)]+[Z(F)]  IV. 1 

forms  the  “complete  impedance  matrix"  of  the  diakopted  antenna,  which 
formulates  the  relationship  between  the  total  terminal  potentials 


M 


produced  by  current  and  field  coupling,  to  all  impressed  currents.  In 
matrix  form 


[I]  =  [Z][I]  IV. 2 

If  the  matrix  elements  1  ( I )  (eq.  III. 19)  and  zM(F)  (eq.  III. 26)  are 
added,  the  resulting  elements  zV1.  have  the  same  formulation  as  those  which 

K  »  J 

pertain  to  field  coupling  between  two  different  elements  (eq.  III. 33).  In 
other  words,  if  the  condition  i  f  i  is  dropped,  equation  III. 33  can  be  used 
as  the  general  formulation  for  all  the  elements  of  the  complete  impendance 
matrix  of  the  diakopted  system. 
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Calculation  of  the  impedances  according  to  eq.  III. 33  requires,  in  principle, 
computation  of  the  scatter  current  and  charge  distributions.  However,  numerical 
results  obtained  with  this  theory  indicate  that  coupling  by  the  scatter  currents 
is  a  negligible  effect.  It  has  been  found  that  coupling  by  the  junction  currents 
prevails  over  field  coupling,  and  field  coupling  by  the  non-Maxwell ian  fields 
dominates  over  that  by  the  (Maxwellian)  scatter  fields.  In  principle,  the  field 
coupling  effect  by  the  scatter  currents  can  be  obtained  with  an  iterative  pro¬ 
cedure  which  is  not  discussed  here. 

If  coupling  by  the  scatter  fields  is  neglected,  the  formula  for  the  elements 
of  the  complete  impedance  matrix  for  the  diakopted  system  reduces  to 


?i»t  _  juj  f  ,xl  'i  i,  ,ci 

Lk,m  v\  Ti  jci^  m*  k  5mqk;d 


k.  r  J  S 
<  m 


1,5.  =  1,2,...,  N 
k  =  1 ,  M . 
m  =  1  ,  M„ 


Thus  all  the  matrix  elements  can  be  computed  from  the  dominant  current  dis¬ 
tributions. 

The  symmetry  of  the  [Z],  i.e. 

zM  =  ZZ’l  IV.  4 

k,m  m,k 

can  be  easily  verified,  by  expressing  in  (III. 33)  the  vector  and  scalar  potentials 
by  the  current  and  charge  distributions  according  to  (III.l)  and  (III. 12). 

Equation  (III. 37)  represents  a  stationary  formulation  of  the  matrix  ele¬ 
ments  of  [Z].  This  means  first  order  errors  in  the  current  and  charge  dis¬ 
tributions  lead  to  second  order  errors  in  the  impedances  (Appendix  4). 
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V.  Interconnection  of  Diakopted  Elements  to  Obtain  Impedance  of  Assembled 

Multi  Element  Antenna 

The  requirement  for  the  diakopted  structure  with  impressed  currents  to  be 
performancewise  identical  with  the  assembled  antenna  are  that: 

a)  The  sum  of  the  impressed  currents  is  zero  at  every  junction  between 
the  structure  elements  and  the  continuity  condition  is  satisfied  at 
every  input  terminal.  This  requirement  assures  that  the  field  of 
assembled  antenna  is  Maxwellian. 

b)  The  scalar  potentials  at  the  interconnected  terminals  are  equal. 

c)  The  potential  difference  between  the  input  in  terminals  is  equated 
with  the  driving  voltage  of  the  antenna  source. 

Imposing  these  junction  conditions;  the  matrix  equation  (IV. 2)  yields  a 
system  of  linear  equations  for  the  unknown  junction  currents  and  the  input  im¬ 
pedance  of  the  antenna.  Using  network  theory  concepts  the  reduction  of  (IV. 2) 
to  this  linear  system  of  equations  by  enforcing  the  junction  conditions  can 
be  formulated  with  a  connection  matrix  [C]  which  reduces  the  number  of  poten¬ 
tials  and  currents  of  the  diakopted  structure  to  those  of  the  actual  structure 
[3].  As  discussed  in  Section  II,  the  impedance  of  the  actual  assembled  antenna 
can  be  written  as: 

[Z]'  =  [C]t[Z][C] 

where  [Z]  and  [Z] *  refers  to  the  actual  and  diakopted  structure  respectively. 


The  following  example  shows  how  [C]  and  [Z] ’  are  obtained. 


Example 


As  an  example  we  apply  the  diakoptic  theory  to  an  ordinary  thin-wire 
dipole  antenna  and  compare  the  results  with  the  exact  data  available  in  the 
literature.  To  obtain  a  multielement  structure  we  cut  each  wire  in  halves,  as 
shown  in  Figure  3,  and  consider  each  half  as  a  structure  element.  The  aiakooteo 
dipole  comprises  two  structure  elements  with  one  terminal  and  two  structure 
elements  with  two  terminals,  so  that  the  total  number  of  terminals  is  six.  The 
complete  impedance  matrix  of  the  diakopted  structure  [Z]  is  therefore  a  5x6 
matrix.  However  there  are  only  3  different  impeoances  because  the  four  struc- 
-.jre  elements  have  been  assumed  to  be  alike. 

using  the  enumerations  of  Figure  2  the  matri  x  equation.  .  1.  .2)  nas  toe  fcrr 


4  =  ^2 


l, 

0 


32  -,23  _  -14  741 

n  ~  ~  Lzz~  Lzz 


7 

“~7 


,43 

‘-21 


If  coupling  by  the  scatter  currents  is  neglected,  the  darkened  pert 
of  the  impedance  matrix  Z  is  the  current  coupling  matrix  Z(I). 

The  interconnection  conditions  require 


where  Iq,  Vq  are  input  current  and  driving  voltage  of  the  antenna. 
Because  of  the  symmetry  of  the  antenna 


Current  and  voltage  matrix  of  the  interconnected  antenna  are 


ll 

fvT! 

[I]' » 

0 

,[♦]’  = 

0  1 

Thus,  the  interconnection  matrix  becomes 
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and  the  impedance  matrix  of  the  assembled  antenna 


For  the  numerical  calculation  of  the  impedances,  the  following  simplifying 
assumptions  have  been  made: 

a)  coupling  by  the  scatter  currents  is  negligible 

b)  the  dominant  current  distributions  which,  in  this  example,  are  the 
same  for  all  the  elements,  can  be  approximated  by  linear  current 
distributions  (uniform  charge  distribution). 

Although  the  latter  approximation  is  rather  crude,  one  should  expect  rea¬ 
sonable  results  if  the  wire  sections  are  short  compared  with  the  wave  length, 
because  all  the  impedance  formulas  are  stationary  expressions.  Linear  current 
distribution  permits  analytic  formulations  of  all  the  impedances  ZQ,  Z-j ,  Z2 
etc.,  and  numerical  calculations  with  a  pocket  calculator  (such  as  HP  25). 

The  results  obtained  are  presented  in  Figure  9.  The  curves  are  plots  (from  a 

table  by  King  [4])  of  the  real  and  the  imaginary  part  of  the  input  impedance 
2L 

of  a  dipole  for  In  —  =  5  as  a  function  of  kL;  2L  is  the  total  length  of  the 

a 

dipole,  and  a  the  wire  radius.  The  crosses  mark  the  values  of  the  input  im¬ 
pedance  from  (45)  with  the  above  assumptions.  For  kL  <  0.8  the  deviation  of 
the  real  part  of  the  input  impedance  from  the  exact  value  is  less  than  10*  and 
for  the  imaginary  part  less  than  1%.  From  this  one  can  conclude  that  the  linear 
approximation  for  the  dominant  current  distribution  is  adequate  if  the  length 
of  a  wire  section  is  <1/1 5a.  This  has  been  born  out  by  computer  results  which 
were  obtained  when  each  dipole  wire  was  diakopted  into  4  equal  sections.  These 
results  are  marked  in  Figure  9  by  dots  and  are  in  good  agreement  with  the  exact 
curves  even  beyond  the  resonance  of  the  antenna. 


VI .  Receiving  Antennas 

In  the  case  of  a  receiving  antenna  the  excitation  is  produced  by  an  ex¬ 
ternal  field  E(e).  The  source  element  is  replaced  by  the  input  impedance  of 
the  receiver.  All  the  quantities  associated  with  the  external  field  should  be 
characterized  by  the  subscript  e. 

The  boundary  conditions  on  the  element  i  yields 


f  r  A  ^ 

jjj«fl(e)  +  70(e)  -  E(e)J*T(ej>dS1  =  o 


i  \ 


Ik  \{e)  =  jw  .  (A(e) - i k  +  <{>(e)qk)dS 

V 


E(e)  T(e)dS1 
M 


where  T^,  q^  are  the  dominant  current  and  charge  distributions  which  the  im¬ 
pressed  current  1^  would  produce  on  this  element.  (e )  the  potential  at  the 
junction  i,k  caused  by  the  external  field.  From  the  boundary  condition  for 
the  dominant  current  distribution  one  obtains 

|  .  +  9rad  3  J'wj  (Ajj(e)  +  S^q (e) )dS*  =  0  VI*2 

c  ^  r*  1 


since  T(e)  is  zero  at  the  junction.  As  shown  in  Appendix  2 


»  ,  ,  ,  *  #  ,  , 
.  (Ak*T(e)  +  ^(eJJdS1  =  . (A(e)*?k  +  $(e)qk)d^ 


Thus  Equation  ( VI . 1 )  reduces  to 


‘k  *k<e)  ■  - {  Wi ds' 
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A  "i 

The  potential  $^(e)  produced  by  an  external  field  at  the  junction  i , k  is  given 
by  the  scalar  product  between  the  external  field  E(e)  and  the  dominant  current 
distribution  function  1 1 »  integrated  over  the  surface  S1  of  the  element.  In 
the  network  presentation  excitation  by  an  external  field  is  equivalent  to  vol¬ 
tage  sources  V(e)  in  series  with  the  terminal  voltages  of  the  intrinsic  imped¬ 
ance  networks. 
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n 


x' 
c 
A 

i  (x) 


p 

U) 


=  wire  radius 

=  ^p,  k  =  I2-  ,  k£  =  B 


VII  Numerical  Results  and  Computer  Programs 
A.  Cylindrical  Wire 

A.l .  Dominant  Current  Distribution,  Dominant  Charge  Distribution  and 
Intrinsic  Impedance  Calculations 
i(x)  =  Current  density 
x  =  Source  point 

=  Observation  point 
=  Velocity  of  light 
=  wave  length 
=  1(o)  =  I0/2*p 
current  distribution  at  x 
l  =  wire  length 


t  « 


/  (x-x1)2  +  4  (£)2  sin2  (a/2) 


The  vector  and  scalar  potentials  at  an  observation  point  x  are: 


1  2u 

Az(*>  -  /  / 

HTr  0  0 


,-jSd 


pdadx' 


$(x) 


J _ 

4ttc£ 


1  2ir 

f  / 

0  0 


1 

-jod 


_d 

dx‘ 


i(x’)  J 


,-JBd 


pdadx1 


Component  of  electric  current  field  intensity  parallel  to  the  surface  of  the 
wiere  is  zero,  and  can  be  written  as: 

*  i  - ° 

or 
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1  2tt  .  j  -j0d  e"j3dn 

;  /  [37.  <(*')  5T'  <1  >  +  si<x'>f  3  d“d>1  =  0 

00 


VII. 1 


sr1(x)  ■  aJcos8(l-xj) 


xj-l<x<xj 


VII. 2 


Xj  =  N  J 


j  =  1,...,  N-l 


Integrating: 


i(x.  ,)  =  i(x.)  +  /  a  cos6(l-x)dx 

J  *  0  y 

a.j 

=  i ( x j )  +  g1  (-sin8(l -Xj_i )  +  Sin6(1 ~Xj ) ) 


1(xN)  =  0 


N  a, 

J  Cg^K-sin6(l-xx_-| )  +  sin8(I-xx))  + 


X=j+1 


-P-  (sine(l-x.)  -  sine(l-x))] 

P  J 

Substituting  VII. 2  and  VII. 3  into  VII. 1  and  choosing 

(i-i) 

x  =  -J-  ,  i  =  1,2,...,N-1,  we  obtain 
N 


VII. 3 


l  a4  (9-ji  +  kii  +  l  ^xji  ^  =  °  1  = 

j=l  J  J1  J1  X=I  AJ 


VII. 4 


j  «j(g-  sin0(l-Xj)  -  j  sin  (l-xJ_1 ))  =  iQ 
J  * 


where 


VII. 5 


9Ji 


,  J  *  d  p-jSd 

[  /  /  cosg(1-x‘  )  (|  )  pdadx‘1 

x.  .  0  ax  a  x'xi 

J  ' 


kji 


[  / 
x 


j-1 


IT 

/  6 

0 


e-jed  , 

d  (sin8(I- 


x.  ,)  -  sinB(l-x))  pdadx1  ]  _  VII. 6 
J  ”  I  *  XX* 


0 

'  Xji 


J  *  -jBd 

[  /  /  3  |  (sin8(1-xx)  -  sin6(l-x>_1))pdadx']  VII. 7 

xj-l  0  x=xj 


Equation  VII. 4  is  simultaneously  solved  to  obtain  a^,...,a^  and  hence  the 
current  and  charge  densities. 

Quantities  VII. 5  to  VII. 7  are  computed  by  quadrature  integration  formula. 
These  expressions  can  be  considerably  simplified  when  ij*j,  resulting  in  com¬ 
putation  saving. 


A. 2  Impedance  Calculations 

Z  =  iiiLL  |=  2vpi  f  i(x’)  A  (x1)  +  i^-,i(x')  f(x')dx'  VII. 8 

o  x=0  Ao  0  z  u  ax 


•  1/2  1  1  n  -jSd 

/  /  /  CenixjKx')  | 
w  V  0  0  0  0  a 


-jBd 


jS)]  dadxdx' 


or 


N  N 


W 


VII. 8 


where 
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iO  'SOLUTION  OF  THE  INTEGRAL  EQUATION  FOR  THE  CURRENT  IN  A  DIPOLE  OR  NON 

OP  OLE 

20  SUBROUTINE  INTEG ( A , B , C , D » LX , MY , X 1 , BETA , S f FUNCTI f RHO , DELTA , THETA ) 

30  REALMS  Zf WEIGHT ,  XI .BETA. XI , YJ, RHO, DELTA, THETA 
40  COMPLEX* 16  FUNC  , S  *  FUNCTI 
50  DIMENSION  Z ( 24 >, WEIGHT ( 24 ) 

60  DATA  Z/0. 577350269 , 0.0,0. 774596669 ,X 

70  0 . 3399S1044  , 0 . 86 1136312. 0,0,0. 5384693 1 0  , X 

SO  0. 906179846  r  0.23861 9186 , 0.661209387,0.932469514,.% 

90  0. 148874339 , 0.433395394, 0.679409568 , 0.865063367,% 

100  0.973906529,  0.0 ,0. 201 194094 ,0. 3941 5: 1347 ,X 

110  0.570972173.0.724417731 , 0.848206583  ,  0.937273392  ,X 

120  0.987992513/ 

130  DATA  WEIGHT/1. 0.0. 388888889. 0.555555556.% 

140  0.652145155.0. 34 7854845.0.568888889. 0.478628671.X 

150  0.23 A 926885.0. 467913935.0, 36 0761573. 0.171324493.0. 295524225 . 0 . 269266 
719.X 

160  0.219086363.0.149451349.X 

170  0.066671344.0.1012891,0,198431485.0.186161000.x 
180  0. 166269206,0. 139570678,0. 107159221 ,0,070366047,X 
L 90  0.030753242/ 

200  S=<0.0,0.0) 

210  DO  10  1=17,24 
220  DO  10  11=1, LX 
230  DO  10  12=1,2 
240  DO  10  J=1 7 , 24 
250  DO  10  Jl=l ,MY 
260  DO  10  J2= 1,2 
270  30  STEPY= ( D-C ) /MY 
230  D1=C+STEPY*J1 
.90  C.l  =  Dl -STEPY 
300  STEPX=(B-A)/LX 
710  B1=A+STEPX*I1 
520  A 1 =B 1 -STEPX 

330  XI* < <-l)*#I2*Z<I)*<Bl-Al)+Bl+Al>/2 
340  YJ=  < (-l>**J2*Z(J)*<Dl-Cl>+ni+Cl>/2 
350  FUNC=FUNCTI (XI, XI »YJ, BETA, RHO, DELTA, THETA) 

360  10  S=S+(B1-A1)*(D1-C1) /4*WEIGHT ( I ) * WEIGHT ( J ) #FUNC 
370  RETURN 
330  END 

390  SUBROUTINE  SIMQ  (A.B.N.KS) 

400  REAL*8  A »  B , BIGA » TOL , SAVE 
410  DIMENSION  A(NtN) ,B(N) 

420  'FORWARD  SOLUTION' 

430  T0L=0 . 0 
440  KS=0 
430  DO  65  J= 1 ,  N 
460  J  Y  =  J+ 1 
470  B I GA=0 
480  DO  30  I=J»N 

>*0  'SEARCH  FOR  MAXIMUM  COEFFICIENT  IN  COLUMN' 

370  IF ( DABS  <  BIGA ) . GE , DABS  <A(I,J)>) GO  TO  30 
IIP  B IGA=A  < I . J  > _ 
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'I 


520  I MAX- I 
530  30  CONTINUE 

540  ’TEST  FOR  PIVOT  LESS  THAN  TOLERANCE  SINGULAR  MATRIX" 

550  IF ( DABS ( BIGA ) >TOL ) GO  TO  40 
560  KS= 1 
570  RETURN 

530  "INTERCHANGE  ROUS  IF  NECESSARY" 

5<?0  40  Ii  =  J+N*(  J-2) 

600  DO  50  K=J,N 
610  SAVE=A  <  J  *  K ) 
o20  A  <  J , K) =A ( I MAX  ,  K ) 

*30  A( IMAX , K ) =SAVE 

•40  "DIVIDE  EQUATION  BY  LEADING  COEFFICIENT" 

<■>50  50  A< J,N)=A< J,K)/BIGA 
660  SAVE  =  B < IMAX ) 

-"O  B<IMAX>=B<J) 

.-,30  B  ( J )  =:SAVE/BIGA 
•  -  0  "ELIMINATE  NEXT  VARIABLE" 

■00  IF  < J=N ) GO  TO  70 
<10  DO  65  IX= JY , N 
7 20  DO  60  JX= JY  » N 

7 30  60  A(IX,JX)=A(IX,JX)-A<IX,J)#A<J,JX> 

740  65  B<IX)=B<IX)-(B< J)*A< IX, J) > 

750  "BACK  SOLUTION" 

760  70  NY=N- 1 

770  IT=N*N 

780  DO  80  J=1,NY 

790  IB=N-J 

800  ION 

810  DO  80  K  =  1 ,  J 

820  B:iB)=B<IB)-A<IB,IC)*B<IC) 

330  30  IC=IC-1 
840  RETURN 
350  END 

360  SUBROUTINE  INTEGR ( INITI A , FINAL , NUMINT , KERNEL , RESUL , XI , X2 , BETA , RHO , DE¬ 
LTA,  THETA) 

870  REAL*8  Z < 7 ) , WEIGHT < 7 ) , INITI A , FINAL , INTERV , A , B , BETA , XI , RHO , X2 , DELTA , T 


HETA 

<380  C0MPLEX*16  S * RESUL , KERNEL 

390  DATA  Z  /O . 201 194094 , 0 . 394 151 347 , % 

900  0.570972173,0.724417731,0. 848206583 ,0. 937273392,% 

910  0.987992518/ 

920  * 

930  " 


940  DATA  WEIGHT  /O. 198431485,0 . 186161000,% 

950  0.166269206,0.139570678,0.107159221,0.070366047,% 

960  0.030753242/ 

970  " 

980  INTERV=< FINAL- INITI A) /NUMINT 
990  A=INITIA 


1 000  RESUL- ( 0 . ODOO , 0 . ODOO ) 
1010  DO  20  J=1 , NUMINT 
1020  B-A+ INTERV 
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1030  ' 

1 o  40  S— < 0 . 0000 1 0 . 0D 00 ) 

1050  DO  10  1=1 t 7 

1. 060  10  S=S+ WEIGHT  <  I )  *  <  KERNEL  ( X 1 » X2  ,  ( Z  ( I  >  *  ( B-A  )  +  B+A ) /2 »  BETA »  RHO  .  DE 

LTA»  THETA ) +KERNEL  <  XI »  X2  ,  <-Z<I)*(B-A)+B+A)/2*  BETA » RHO »  DELTA » THETA  >  ) 

1070  S=S+0.202578242*KERNEL(X1»X2  ,<B+A>/2>  BETA . RHO > DELTA » THETA > 

1080  S=(B-A)/2.0*S 

1090  A=A+INTERO 

1100  20  RE3UL=RESUL+S 

LI  10  RETURN 

1120  END 

1130  • 

1140  * 

1150  FUNCTION  G < XI > X2 » ALPHA » BETA » RHO , DELTA , THETA ) 

1160  COMPLEX* 16  GrE,El 

1170  REAL*8  XI »  X2*  BETA* RHO » ALPHA  tDrDl * DELTA*  THETA 
1.180  D=DSQRT  (  <  XI -X2>**2+<  2. 0*RH0*DSIN<  ALPHA/2.0)  >**2> 

1190  D1  =  DSQRT < < XI +X2+DELTA ) **2  + ( 2 . 0*RH0*DS IN ( ALPHA/2 . 0 ) )**2> 

1200  E=  < 1 . ODOO  »  0 . ODOO ) *DCOS ( BET A*D ) - <  0 . ODOO  1 1 . ODOO ) *DSIN ( BETA*D ) 

1210  El=( 1 , ODOO *0. ODOO )*DC0S<BETA*D1 >-<0.0D00* 1 . ODOO > *DSIN ( BETA*D1 ) 

.1 220  G- (  -E/D**2-  <  0 .  ODOO  r  1 .  ODOO  >  #BETA*E/D )  *  ( XI ~X2  >  /D-THETA*  (  -El  /D 1  *#2-  (  0  . 

ODOO f 1 . ODOO ) *BETA*E 1/DI ) * < X 1+X2 ) /D1 

1230  RETURN 

1240  END 

1250  * 

1260  FUNCTION  H(X1 * X2 *  ALPHA  *  BETA  * RHQ *  DELTA *  THETA > 

1270  COMPLEX* 16  H*E,E1 

1230  REAL*8  X 1 * X2 * BETA » RHO , ALPHA * D > Dl * DELTA * THETA 
1290  D=DSQRT< (Xl-X2)**2+(2.0*RH0*DSIN<ALF'HA/2.0) )**2) 

.1  300  Dl-DSQRT  (  ( X 1 +X2+DELTA )  **2+  <  2 . 0*RH0*DSIN  ( ALPHA/2 . 0 )  )  **2 ) 

1310  E= ( 1 . ODOO * 0 . ODOO ) #DCOS  <  BETA*D ) - <  0 . ODOO  » 1 . ODOO ) *DS I N ( BETA*D ) 

1320  El  =  < 1 .0D00,0. ODOO )*DCOS< BETA*01 >-<0.0D00* 1 . ODOO) *DSIN ( BETA*D1 ) 

1330  H=E/D+THETA*E1/D1 
1340  RETURN 
1350  END 

1551  FUNCTION  HPR ( ALPHA * XI * X2 * BET A * RHO * DELTA r THETA ) 

1:52  COMPLEX* 1 6  HPR,E»E1 

1353  REAL*8  XI > X2 * BETA * RHO * ALPHA * D * D 1 * DELTA * THETA 

1354  0=DSGRT< ( XI -X2 > **2+ ( 2 . 0*RH0*DS IN < ALPHA/2 . 0 > )**2> 

1355  D1=DSQRT ( ( X 1 +X2+DELTA ) **2+ < 2 . 0*RH0*DS IN < ALPHA/2 . 0 ) >**2> 

1 356  E= ( 1 . ODOO * 0 . ODOO ) *DCOS ( BETA*0 ) - ( 0 . ODOO » 1 . ODOO ) *DS IN ( BETA*D ) 

1357  El=( 1 . ODOO *0. ODOO >*DC0S(BETA*D1 >-( 0 . ODOO *  1 . ODOO ) *DSIN< BETA*D1 ) 

1358  HPR=E/D+THETA*E1/D1 

1359  RETURN 
l  560  END 

1361  FUNCTION  HI < ALPHA > XI  * X2 *  BETA » RHO » DELTA  *  THETA  ) 

1370  C0MF'LEX*16  HlfE.El 

1380  REAL*S  XI ?X2r BETA »RHO» ALPHA »D*D1 » DELTA*  THETA 
1390  D-DSQRT  < < X1-X2 > **2+ < 2 . 0*RH0*BSIN< ALPHA/2 . 0 ) )**2) 

1400  D1=DSQRT < <Xl+X2+DELTA>**2+< 2, 0*RH0*DSIN< ALPHA/2.0) )**2) 

1410  E=( 1 . ODOO *0. ODOO )*DCOS<BETA*D)-<0. ODOO* 1 . ODOO ) *DS IN < BETA*D > 

1420  El  =  < 1 .ODOO*O.ODOO)*DCOS<BETA*D1 >-<0.0D00* 1 . ODOO ' *DSIN < BFTA*Dl ) 

1430  HI =E/D-THETA*El/Dl  +  < 1 . 0-THETA )*< 0 . ODOO *  1 .0000  >*BETA 
L44C  RETURN 


1450  END 

1  4<s0  FUNCT I  ON  IMPKER  (  XI  -  X2  -  ALPHA »  BETA ,  RHQ  -  DELTA  -  THETA  ) 

1470  C0MPLEX*16  IMPKER »E -El 

1480  REALMS  X 1  - X2 -  ALPHA  »  BET  A  -  RHO -  D  -  DELTA  -  THETA  »  D.l 
1490  D=DSQRT ( ( XI -X2 ) **2  + ( 2 , 0*RH0*DS I N ( ALPHA/2 . 0 ) )**2) 

1500  D1=DSQRT  < < X 1  +  X2+DELTA ) **2+ < 2 . 0*RHQ*DSIN ( ALPHA/2 . 0 > >**2) 

1510  E=  ( 1  . ODOO  - 0 . ODOO ) *DCOS ( BETA*D )  -  (  0 . ODOO - 1 . ODOO ) *DSIN  (  BETA*D ) 

1520  El=( 1 . ODOO, 0. ODOO )*DC0S(BETA*D1 >-(0.0D00- 1 . ODOO ) *DSIN( BETA*D1 ) 

1 530  I MPKER=DCGS ( BETA*  ( 1 . 0-X2 ) ) * ( E/D-THETA*E1/D1  +  ( 1 . 0- THETA ) * ( 0 . ODOO -1.0 
DOO)*BETA) 

1540  RETURN 
1350  END 

1560  EXTERNAL  G - H - IMPKER - HI , HPR 

1570  C0MPLEX*16  RESG-RESH-G-H-A<40-40) - AB  <  40  > -  RE SUL  -  CHAR ( 30  > »  CURR - IMPKER 
••  Z  -  SUMH ,  ZPR  -  RESU1  -  CUR  ( 30  > -  IMP , HI » RESH1 - CHARGE »  HPR 

1 530  REAL*8  XI  - X2 ,  STEP  -  BETA , RHO -  X - REALPA , IMAGPA - REAATB - IMAATB - RA ( 1600  >  - R 
Ti  <  50 )  , ATA  ( 1600  )  -  ATB  <  50  )  - MAGCHA  ,  MAGCUR  -  IDCHA  -  IDCUR -  DELTA  ,  THETA  -  ALPHA 
1590  WRITE  (6-2200) 

1600  READ  (5-*)  N- BETA, RHO -DELTA -THETA 
1610  IF  (THETA-O. 5)  171-171,172 
1620  171  WRITE  (6-1171) 

1630  GO  TO  173 

1640  172  WRITE  (6-1172) 

1650  173  WRITE  (6-999)  BETA-RHO-DELTA 

1660  STEP=1 . 0/N 

1670  NTIM2=2*N 

1680  NMINl-N-l 

1690  X1=STEF 

1700  DO  10  1=1 -NMIN1 

1710  X2=STEP/2 . 0 

! 720  SUMH= (0. 0-0.0) 

1730  DO  20  J=1 - N 

1740  CALL  INTEGR (0.0 -3. 1415-1 -G-RESG- XI -X2- BETA, RHO -DELTA- THETA) 

1750  CALL  INTEGR ( 0. 0 - 3. 1 415 - 1 -H - RESH- XI -X2- BETA -RHO -DELTA -THETA) 

1760  SUMH=SUMH+RESH 

1770  A(I-J)  =  (DSIN( BET  A* ( 1 » 0— ( X2-STEP/2 . 0 ) ) ) -DSIN ( BETA* ( 1 . 0- ( X2+STEP/2 . 0 ) 

)) )/BETA*RESG+( DSIN ( BETA* (1.0-(X2~STEP/2.0> ) >-DSIN(BETA*(1.0-X2> ) )*BETA/ 

N*RESH+ ( DS IN ( BET  A* ( 1 ♦ 0- ( X2+STEP/2 . 0 ) ) ) -DSIN ( BETA* ( 1 . 0-( X2-STEP/2 .0  > ) ) )*B 

ETA/N*SUMH 

1780  20  X2=X2+STEP 

1790  10  X1=X1+STEP 

1800  X=STEP 

1310  DO  30  J=1  -N 

1820  A(N,J)=( DSIN ( BETA* ( 1 . 0-X ) ) -DSIN ( BETA* ( 1 . 0- ( X-STEP ) > ) >/BETA 

1330  30  X=X+STEP 

1S40  DO  40  1=1, N 

1350  40  AB( I )  =  I/N*1 ,0 

1860  DO  50  1=1, N 

1870  DO  50  J  =  1  -  N 

1380  REALPA= ( A ( I  - J ) +DCON JG ( A ( I »  J ) )>/2.0D00 

1390  IMAGF'A=  ( A  ( I  -  J  )  -DCON  JG  <  A  ( I  -  J )  )  ) /2 .  ODOO/ ( 0 .  ODOO  - 1 .  ODOO  ) 

1900  REAATB= ( AB ( I )+DCONJG(AB( I ) ) )/2.0D00 

i'UO  IMAATB  =  ( AB  ( I )  -DCONJG  ( AB  ( I )  )  ) /2 .  ODOO/ (  0  .  ODOO  - 1  .  ODOO  ) 

1920  RA( I+< J-l )*NTIM2)=REALPA 


L >30  RA ( I  +  C J+N-l )*NTIM2)=-IMAGPA 

1940  RA< I+N+C J-l )*NTIM2)=IMAGPA 

1950  RAC I+N+C J+N-l ) *NT IN2 > =REALPA 

1960  RBC  I )=REAATB 

1970  RBC I+N)=IMAATB 

1980  50  CONTINUE 

1990  DO  410  1=1 fNTIM2 

7000  DO  410  J=1 , NTIM2 

2 010  ATAC ( J-l )*NTIM2+I )=0.0D00 

2020  DO  410  K=1fNTIM2 

2030  410  ATAC ( J-l ) *NTIM2+I ) =ATA ( C J-l > *NT IM2+ I > +RA C ( I - 1 ) *N*2+K ) *RA ( ( J-l )* 
N*2+K) 

2040  DO  420  1=1 f NTIM2 
2050  ATB ( I ) =0  . ODOO 
2060  DO  420  K=1fNTIM2 

2070  420  ATB  C I ) =ATB CI)+RACCI-1)  *N*2+K ) *RB  C  K ) 

2080  DO  430  I=1fNTIM2 
2090  430  RB < I ) =ATB Cl) 

2100  CALL  SIMC)CATAfRBfNTIM2fKS) 

2110  WRITE  ( 6 f 3000 ) 

2120  WRITE  C  6  f 100 ) 

2130  DO  110  1=1 fN 

2140  110  WRITE  <6f*)  RBC I ) fRBC I+N) 

2150  WRITE  <6,100) 

2160  WRITE  C  6  f  2000 ) 

2165  WRITE  C 6  f 6000 ) 

2170  X=ST EP/2.0 
2180  DO  150  I=1fN 

2190  CHAR C I ) = C RB C I ) * C 1 ♦ ODOO f 0 ♦ ODOO ) +RB C I+N ) # C 0 . ODOO f l . ODOO ) > *DCOS C BET A* C 
1 .O-X) ) 

2195  CHARGE=CHAR  C I ) /BETA 
2200  MAGCHA=CDABS  C  CHAR ( I )  ) 

2210  MAGCHA=MAGCHA/BETA 

2220  WRITE  C6f5000)  X f CHARGE f MAGCHA 

2230  150  X=X+STEP 

2240  WRITE  C6f100> 

2250  WRITE  <6f2100) 

2255  WRITE  (6,6000) 

2260  CURR= CO. ODOO fO. ODOO) 

2270  X=1 .0 

2280  DO  160  1=1 fN 

2290  CURR=CURR+  C DSIN C BETA# C 1 . O-X ) ) -DSIN C BETA* C 1 . 0- C X-STEP ) ) ) ) /BETA* C RB C N 
-I  + 1 ) *  C 1 . ODOO  f  0 . ODOO ) +RB C  N- 1  + 1 +N ) *  C  0 . ODOO  f 1 . ODOO  )  ) 

2300  X=X-STEP 

2310  MAGCUR=CDABS  C  CURR ) 

2320  CUR C N-I  +  l >  =CURR 

2330  160  'WRITE  <6f5000)  X  f  CURR  f  MAGCUR 
2340  WRITE  <6f100> 

2350  WRITE  C6f2300) 

2360  *  ZPR  IS  THE  POTENTIAL  AT  THE  END 
2370  Z= C 0. ODOO fO. ODOO) 

2380  ZPR= CO. ODOO fO. ODOO) 

2390  X2=0 . ODOO 
2400  DC)  170  1  =  1.  f N 


47 


2-410  CALL  INTEG  (  X2  » X2  +  STEP  f0.0»3.1415*l»l»0.0»  BETA » RESUL  » I MPKER »  RHO ,  DEL  T 
A. THETA) 

2  *20  CALL  INTEG  (  X2 » X2+STEP  * 0 . 0r3.1415tlrlrl.0r  BETAr  RESU1 r I  MPKER  r  RHO  r  BELT 
h  *  THETA ) 

2  430  Z=Z+RESUL* <  RB  < I >  *  < 1 . ODOO r 0 . ODOO ) +RB < I +N ) *  <  0 . ODOO  r 1 . ODOO ) ) 

2440  ZPR=ZPR+RESU1 * < RB < I ) * ( 1 . ODOO r 0 . ODOO  >  +RB <I+N)*(0. ODOO r 1 . ODOO ) ) 

2450  170  X2=X2+ST£P 

2430  Z=Z*377,0#<0. ODOO, 1. ODOO )/4. 0/3. 1415**2/BETA 
2470  ZPR=ZPR*377. 0*<0. ODOO r  1. ODOO )/4. 0/3. 1415**2/BETA 
2480  ALPHA=1 . 0471 
2490  IMP= (0. ODOO >0. ODOO) 

2500  Xl=STEP/2»0 
2510  DO  66  J=lrN 
2520  X2=STEP/2 . 0 
2530  DO  55  1=1  r  N 

2540  D=DSQRT( < XI -X2 ) **2+ < 2 . 0#RH0*DSIN ( ALPHA/2 . 0 ) )**2> 

2550  D1=DSQRT< < X1+X2+DELTA ) **2  + ( 2 . 0*RHQ*DSIN< ALPHA/2 . 0 ) )**2> 

2560  E=( 1 .0D00r0.0D00>*DC0S<BETA*D)-(0.0D00r 1 . ODOO ) *DSIN < BETA*D > 

2570  El  =  <  1  . ODOO r 0 . ODOO ) *DCOS  <  BETA*D1 ) - <  0 . ODOO  r 1 . ODOO ) *DSIN < BETA*D1 ) 

2572  CALL  INTEG ( Xl-STEP/2 . 0 r X 1 +STEP/2 . 0 r X2-STEP/2 . 0 r X2+STEP/2 . 0 r 1  *  1 » 0.78 
54 i BETA  f RESH  r  HPR r RHO  r DELTA r  THETA  > 

2573  CALL  INTEG  <  Xl-STEP/2. 0  r  X 1 +STEP/2 . 0  r  X2-STEP/2 . 0 r  X2+STEP/2 .  0  r 1 r 1 r  0 . 78 
54 ,  BETA r RESH  1  r  HI  r RHO ,  DELTA r  THETA ) 

2580  IMP=  IMP+BETA**2*CUR  <  I )  *CUR  ( J )  ♦RESH-  (  ( 1 .  ODOO  r  0 .  ODOO )  *RB  ( I )  +  (  0  .  ODOO  r 1 
. ODOO ) *RB ( I  +  N ) ) * <  < 1 . ODOO r 0 . ODOO ) *RB <  J )  +  ( 0 . ODOO r 1 . ODOO ) *RB ( J+N ) ) *DCOS ( BET 
A* < 1 . 0D00-X1 ) ) *DCOS ( BETA*  < 1 . 0D00-X2  > ) *RESH1 
2590  55  X2=X2+STEP 
2600  66  X1=X1+STEP 

2610  IMP=IMP*< 0. ODOO, 1. ODOO) #377.0/4. 0/3. 1415/BETA* ( 1 .ODOO+THETA) 

2640  WRITE  <6r*)  IMP 
2650  WRITE  <6r2400) 

2660  WRITE  <6r*>  ZPR 
2670  100  FORMAT ( '  ') 

2630  2000  FORMAT ( '  DISTANCE  CHARGE') 

2890  2100  FORMAT ( '  DISTANCE  CURRENT') 

2693  6000  FORMAT ( '  DISTANCE  REAL  IMAG  MAGN ' ) 

2700  2200  FORMAT ('  N r BETA r RHO r DELTA r THETA ' ) 

2710  2300  FORMAT ( '  IMPEDANCE') 

2720  5000  FORMAT ( F6 . 3 » 3F8 . 3 ) 

2730  3000  FORMAT ('  ') 

2740  2400  FORMAT ( '  POTENTIAL  AT  THE  END') 

2750  1171  FORMAT ('  MONOFOLE') 

2760  1172  FORMAT <'  DIPOLE') 

2770  999  FORMAT ( '  BETA= ' r F8 . 4 r ' RADIUS* ' r F8 . 4 r ' GAP= ' r F8 . 4 > 

2 “30  STOP 
2790  END 
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VIII.  Dominant  Current  Distribution  and  Impedance  of  a  Circular  Disc  Fed 


Integrating 


N 


KxJ  =  l  a,  (/TF  -  /i-xf  J 

K  A=|<+1  A  A  A- 1 


k  =  0,  ....  M-l 


I(X)  =  I(xk)  +  ak+1  (/l-x*  -  /1-X2)  xk<'x$xk+1 


Equation  VI 1 1. 4  can  be  reworked  as: 


N  N 

l  l>.(g-(x)  +  k,(x))  +  l  a.£  ,(x)]  =  0  0^x<1 

1  =  1  J  J  J  ^=j  A  AJ 


j=l 


where 


x . 


j  *  d  ,e“jed 


9j(x)  •  /  /  fj(f  )(-==)  dadx1 

J  x .  i  0  ax  a 

J  “  ■ 


X  . 


kj(x) 


j  -jsd  _ _ 

/  J  (32  rr  cosa  (/1-X?  , 
v  n  u  J* 1 

xj-i  u 


-  A-x2  )  dadx‘ 


X-  TT  -j6d 
■  J  f  a2  e 


*Aj(x)  =  /  J  /  s2  ^  cosa(/l-x2  -  /l-x2_i )  dadx* 


x.  ,  0 


Rewriting  VI 1 1.5  and  adding  the  terminal  condition, 
N  j 


jl,  aj  (90(xi>  *  k0<Xi)  *  J,  W’ 


=  0  x,  =  1-1/2 


i  »  1,  ....  N 


jl,  aj  -  /rTTT)  ■  '< 
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#*»«*•**"  '  - 


-  i*  •  i Cl#  - 1  ■■  «  * 


VIII. 5 


■fW 


Solution  of  above  equations  yields  the  dominant  current  distribution. 
Impedance  is  computed  as 


7  -  $(x' )  I 
L  ~  T  I 

Lo  x‘=0 


-  N  x, 

_ J —  f  a  r  k 

4rreRai  .  f  ak  J 
k-i  xk_1 
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Program  for  Current  and  Impedance  Calculation 
of  a  Center  Fed  Circular  Disk 


y  •.  }  •  ;“•»  A  '  i  i  „  ...  .  .  :>  J  .  >'■  H  v  J  /  w  X 

*.  %  i A  J  *v  av-  %  J  .■  f  r  i  \  .  \  ; 


,-XS  P-A&5  LS  BSST  ftUALlXI  ' 

■*m  oki n  jischlshsd  so  8M  - - .* 


T 


IX.  Impedance  Calculation  of  a  Thin  Wire  with  Linear  Current  Distribution 

Linear  Current  distribution 

=5 


,-M  1  1  *  -jBd  A  ,  -  j3d 

z  =  £  J  /  C32i(x)i(x*)  |  -fc(x)fr(x')(|  +  j$)  dxdx'da 

M.ij/i 

4uz  e 

Assumption  of  linear  current  distribution  and  computation  of  the  kernel  at 
o  *  n/3 

iM  1  1  e'J*6d  e‘J’gd 

1  =  S  I  I  [B2(l-x)(l-x')  |  -  (|  +  36)]  dxdx‘ 


0  0 


d  =  /  (x-x‘ )z  +  p2 
p  =  Radius  of  the  wire 
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Program  for  Impedance 


of  a  Wire  with  Linear  Current  Distribution 
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} 


'&.US 


**'■*  LS  Hljf  %U4i>lTv  i.j.A  .. 

■-*»£*  naeuLsjgp  » 


»PJ*S 
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Program  to  Compute  Impedance  of  Dipole  Assuming  Linear  Current  Distribution 


I'rfis  Pia*  u  Basr  mum  aymmja* 
oon  nwiisfl®  m 
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I  » 
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pa&A*i^ 


XI. 


Computation  of  Dominant  Current  Distribution  for  all  Frequencies  Via 
Static  Charge  Distribution 

Accurate  computation  of  dominant  current  becomes  one  of  the  most  important 
tasks  in  using  Diakoptic  Theory  for  complicated  radiating  structure  analysis. 
For  a  general  structure  element  it  is  a  difficult  task.  Furthermore,  if  we  can 
compute  dominant  current  for  all  frequencies,  the  impedance  spectrum  becomes 
easy  to  compute.  In  many  shapes,  such  as  spheres,  circular  disks  and  cylindri¬ 
cal  conductors,  the  static  charge  distribution  is  either  known  or  easy  to  com¬ 
pute.  In  this  section,  we  shall  develop  an  algorithm  to  compute  dominant 
current  distribution  from  static  charge  distribution. 

Let 

T(r')  =  l  (-jk)n  Or' )  ,  i  (r* )  ■  T'  XI. 1 

0  n  n 

q(r’)  *  l  (-Jk)U  qn(r‘)  ,  qn(r')  =  qp  XI. 2 


q '  =  V  -  i  * 

^n  n  ' 


k  = 


Thus 


00  GO 


*  (?)  ■  £  l  I  <-jk)n+"  *nX(r) 

4  n=0  m=0  nA 


4tt 


nX 

Similarly 


1  \  (-jk)x  I, 

X=0  n=0 

T  (7* )Dx-n_1 
tji  ( X-n ) 


f  in(r‘)D; 

i,  (X-n)! 


hX 


dS*  ,  D  = 


r  -  r 


00  X 

+<?>  -  ET  l.  l. 


X=0  n=0 


where 


XI. 3 


XI. 4 
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V<r>  ■  !  Tx^nfr  °x'n''<iS' 

S 

S'  represents  the  total  structure  area. 


Taking  gradient  of  XI. 4 

oo  ^ 

-jwye  =  J-  l  l  (-jk)A  T  .  (r) 
4ir  X=0  n=0  nA 

where 


^  °A'n‘z  If  <  ds' 

S 

Boundary  Condition  Etan  =  0,  implies 

(-k2A(7)  +  ju)ueV(£)xd'5  =0,  d^  =  rfdS 


Substituting  XI. 3,  XI. 4  into  XI. 5  and  equating  powers  of  (-jk). 


0  .  /  mr'i'  3d  rxn  ,c ,  _  n 

*00  % !  (V*1o)p-3FH'dS  =0 


p  =  [  f y.i '  )  _L  Id  IHL  He '  =  0 
*11  V  Dz  3r  ab  u 


f(V*1X+l)  W  3?  ™  dS'  =  JQ  ('*n,X+l  +  An,X-lXn)  X=1,2’"-  XI*8 

n  £  X 

If  the  impressed  current  is  assumed  real  at  all  frequencies,  we  obtain 
i^  =  0,  ik(o)  *  0,  k  =  2,  ...,  n.  Equations  XI. 6  and  XI. 8  compute  the 
static  current  iQ  and  rest  of  the  currents  i^»  ij>  ...»  in. 


A.  Cylindrical  Conductor 

Dominant  current  equations  are  simplified  as: 

\  \  (V*Tq(x‘  ))  ^^-dx'dct  =  0 

0  0  0  D3 


2irpi0(o)  =  IQ  ,  D  =  /(x-x1  )2  +  ^psinct/^)* 
i7(x)  =  0 


v  1 

/  /  (v-T 
0  0 


•'  J  dx'da 

x+r  D3 


X-l  tt  1  (V*i’'(X-n)(x-x')-i'  ,  ? 

I  1  1  c - JpsTTT! - V"2dx'dc, 

n=0  0  0  1  n  U- 


Tables  1  and  2  represent  computed  values  of  the  current  In( x)  and  charge  qn(x). 
Figures  19  an<^'20  show  current  In( x)  and  charge  qn ( x)  for  different  values  of  6 


9 

i„(x> 

Table  1. 

of  a  cylindrical 

21 

conductor  Q  =  2$,n  —  = 
P 

i'i 

10  • ' 

I0(x) 

I2(x) 

I3(x) 

Vx) 

X 

W 

w 

y®) 

►— « 
O 

O 

I  0.0 

1.000 

0.0 

0.0 

0.0 

1  0.1 

0.865 

-0.0644 

-0.0070 

0.0095 

1  0-2 

0.768 

-0.1000 

-0.0110 

0.0157 

I  0.3 

0.677 

-0.1231 

-0.01 38 

0.0206 

!  o-4 

{ 

0.588 

-0.1357 

-0.0155 

0.0243  1 

0.5 

0.500 

-0.1386 

-0.0163 

0.0264  j] 

0.6 

0.412 

-0.1323 

-0.0160 

0.0266  U 

1  •' 

0.7 

0.323 

-0.1169 

-0.0146 

0.0247 

•  3 

0.8 

0.232 

-0.0923 

-0.0120 

0.0203  [ 

0.9 

0.135 

-0.0570 

-0.0079 

0.0141  U 

1  8 

T  .0 

0.0 

0.0 

0.0 

0.0  jj 

; 

i 

i 

:  ; 

' 

! 

t 

i 

; 

« ^ 

- 

I 
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Table  2.  :  y-^y  •  — ^ —  of  a  cylindrical  conductor  £1  =  2£n  y-  =  10 


X 

1  «o<*> 

I0(o)  dx 

■j  dl2(x) 

ysr  dx 

,  dl-(x) 

yol*  dx 

1  dl^(x) 

yoy*  dx 

0.1 

-1.010 

-0.427 

-0.0466 

0.0683 

0.2 

-0.929 

-0.286 

-0.0331 

0.0557 

0.3 

-0.896 

-0.176 

-0.0225 

0.0435 

0.4 

-0.881 

-0.076 

-0.0125 

0.0292 

0.5 

-0.877 

0.018 

-0.0024 

0.0124 

0.6 

-0.881 

0.108 

-0.0081 

-0.0073 

0.7 

-0.896 

0.199 

-0.0195 

-0.0304 

0.8 

-0.929 

0.293 

-0.0325 

-0.0579 

0.9 

-1.010 

0.414 

-0.0500 

-0.0957 
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x  =  r/Jt 


1.0 


FIGURE  19  Charge  Distribution  on  a  Cylindrical  Conductor 
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B.  Circular  Disc 
2irxi(x)  =  I(x) 

D2  =  (x2  +  x'2  -  2xx'cosa)1//2 
x  =  r/R,  x'  =  r'/R 


Static  charge  equation  is 

<*'';>  ST  <ff>  dx'da  a  0 


2TTXi(x)  =  I0 


i^x)  =  0 

Higher  order  current  densities  are  obtained  by 


i’  x*  .  9 

-  >  cosa]  Da  d  dx'da 

Lim  i  (x)  -  0  X  =  1,2,... 

x-K)  A 


Tables  3  and  4  show  computed  values  for  various  currents  and  charges. 
Figures  21  and  22  show  current  distribution  and  charges  for  different  values 
of  3. 
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Exact 


Exact 

or*0 

Static 

/1-x2 

Vx> 

I2(x) 

I3(x) 

I4(x) 

w 

yur 

ysr 

1.000 

1.000 

0.0 

0.0 

0.0 

0.995 

0.995 

-0.014 

-0.001 

0.00030 

0.980 

0.980 

-0.039 

-0.006 

0.00105 

0.954 

0.954 

-0.070 

-0.012 

0.00188 

1 

0.917 

0.017 

-0.101 

-0.021 

0. 00241 

0.866 

0.866 

-0.129 

-0.031 

0.00229 

0.800 

0.801 

-0.150 

-0.041 

0.00125 

0.714 

0.715 

-0.161 

-0.049 

-0.00074 

0.600 

0.600 

-0.157 

-0.054 

-0.00333 

0.436 

0.433 

-0.129 

-0.049 

-0.00526 

0.0 

0.0 

0.0 

0.0 

0.0 

Table  3  : 


v*> 


of  a  circular  plate  fed  at  the  center 
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Exact 


X 

Static 

-1 

2Wl-x2 

i0(°) 

I„(») 

v°> 

v°> 

0.0 

-0.159 

-0.159 

-0.647 

-0.045 

+0.011 

0.1 

-0.160 

-0.160 

-0.330 

-0.045 

00.009 

0.2 

-0.162 

i 

-0.162 

-0.229 

-0. 043 

!  0. 007 

: 

0.3 

-0.167 

-0.166 

-0.165 

-0.041 

0. 004 

0.4  1 

1 

!  -0.174 

1 

-0.173 

-0.117 

-0.037 

0.001 

0.5 

-0.184 

-0.182 

-0.081 

-0.033 

-0.002 

0.6 

-0.199 

-0.198 

-0.043 

-0.026 

-0.004 

0.7 

-0.223 

-0.222 

-0.011 

-0.017 

-0.006 

0.8 

-0.265 

-0.262 

0.025 

-0.004 

-0.005 

0.9 

-0.365 

-0.362 

0.082 

0.022 

-0.00185 

Table  4:  d^v  of  a  circular  plate  fed  at  the  center 


div  i ( x ) 


; _ i  _ ; _ _ _ _ _ _ _ _  _  _  ' . . .  _  _ _ 

1  - .  . •  •  -* - * - g  =  1.9 


Figure  21  Charge  distribution  on  a  circular  plate  fed  at  the  center 
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XII.  Top  Loaded  Dipole  Antenna 

In  this  section  we  shall  compute  impedance 
dipole  as  shown  in  Figure  XI-1 

R 


characteristics 


of  a  top  loaded 


y 

Element  3 

3  ^ 

Terminal  {^) 

Terminal  (^) 
Element  4 


T*F\  1 

^■Terminal  (-j) 
2p  * -  Element  1 

_ --Terminal  (1) 

— *-l4  * 

ttk.  2 

Terminal  (^) 
Element  2 
— -Terminal  (^) 


Fig.  22  Dipole  with  Circular  Plates  Capactive  Loading 

Section  X  is  used  to  compute  dominant  current  distribution  on  various  radiating 
elements.  Using  the  dominant  current  distributions,  the  various  impedances  are 
computed  as 


i,  K  ■ 


XII-1 


where  (^)  represents  the  impressed  current  terminal  and  (*)  represents  the 
terminal  where  resulting  potential  is  computed. 
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RUTGERS  -  THE  STATE  UNIV  PISCATAWAY  NJ  COLL  OF  ENGJN-- ETC  F/G  9/5 
MATHEMATICAL  MODELING  OF  MULTI-ELEMENT  MONOPOLE  ANTENNAS. <U> 

MAR  81  N  N  PURlf  G  GOUBAU*  A  STAVRIOIS  DAAG29-79-C-0201 

ARO- 154 15. 2-EL  NL 


■o' 


(z'-2z;tz0+2z6-z6-z7)v 


2(Z„-Z,)(Zi-2Zi*Z0*2Z6-Z5-Z7>-(Z3-Z2-Z4tZ3>(2Z3t2Z3-2Z2-2Z4) 


Thus 


Z.  =  _£ 
in  !o 


_  2(VZ1 ' ■Zo~2Z1+Zo*2Z6~Z5~Z7^~^Z3~Z2~Z4+Z  j^Z3+2Z3~2Z2~2Z4^ 
(Z;-2ZiV2Z6-Z5-Z7>  "  . . ‘ 


I  i'1  .i.  w  i"  ;..J  U  Fs.  HI  I  wJi'ir  L)  I  Li.w  i  he  CiC.L.  I‘  Xili"  Li.liihLL  J  Mw  v^uLL.  i'i  J  i  1  j  L. 
iiUTU  \L  IMPEDANCES  CF  THE  ELEMENT  a  CF  Hi!  ANiEuNii  CuH'Sii  r*NG  Ul: 
A  DIPOLE  WITH  END  CIRCULAR  P  LA  TEC. 

THE  OLTPU  i  OF  THIS  PRO GRAN  AKL  GEL!"  AND  ■  1 U  i  CAL  l.'l  'EDi  nLL  j 
OF  THE  ELEMENTS »  THE  IMPEDANCE  CF  THE  dIPjL-  U1.MC-:  i .  ii_ 
PLATES  a  THE  IMPEDANCE  OF  THE  DIPOLE  Ullll  THE  FLhiES  A  Mi' 


PLATES » THE  IMPEDANCE  OF  THE  DIPOLE  Ull  ll  THE  FLhiLS  AMI' 

THE  CURRENT  AND  CHARGE  DISTRIBUTIONS  ON  I  HE  CYLINDRICAL 
CONDUCTORS  AND  THE  PLATES. 

THE  DOMINANT  CURRENT  AND  CHARGE  DISTRIBUTIONS  OF  THE  CIRCULAR 
PL  ATE  AND  THE  CYLINDRICAL  CONDUCTOR  ARE  LOMPU1 ED  IN  TUG  SEPaRA 1 1 
PROGRAMS  AND  MUST  BE  SUPPLIED  AC  INPUi  DA  iA  TO  THIS  7  RCGR.-.M 
FUNCTION  rJ<XfJ> 

C  CNMON  .ABC  /  CUR VEC  <  5  /  1 1  >  y  CHAVEC  Hi  y  I  i  -  >  CURVED  A»iut  Cl  IAVEd  >j»U) 
•COMMON  /DEF/  F.HC  ■*  R  t  BET  A 
COMMON  /GHI/  I»M,K 

!  .■  M  i-  I.  L  \  1  w  i"  W  y  C  <J  !\  V  E  w  y  Ci  l  A  V  E  L.  y  C  L  K  V  u. *.*  y  u  i  !  r  i  V  u.  D 
L  -  10,04tXfi 
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10  CONTINUE 
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RETURN 

END 

FUNCTION  FIKXy  J) 
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*  G  i  t  »-i  i./(0.  \,  y  1  . 
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f  i  r*  r  r  rr 
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(6/5103)  RES 10 

FORMA 

T( '  24  «  23122  a 

'  /2F.L4 

*  vj  .' 

ZDW  r2 

.  ODOQ-K  (RES-RES1  .• 
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(3 y 3i:o;  ZDW 
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l  /I  r  E  1 1 1  i  i  2  l.  Ui 
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123  CONTINUE 
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Conclusion 

The  major  advantage  of  the  diakoptic  theory  for  multielement  antennas, 

is  that  the  problem  of  determining  the  current 
distribution  on  the  antenna  need  not  be  solved  for  the  structure  as  a  whole, 
but  only  for  the  Individual  structure  elements.  Excitation  of  each  structure 
element  is  ascribed  to  the  currents  at  its  junction  with  adjacent  elements, 
and  to  the  fields  of  the  surface  currents  on  all  the  other  elements.  The 
current  distributions  produced  by  the  junction  currents  have  been  termed 
dominant  current  distributions,  because  they  constitute  the  major  portion 
of  the  currents  on  the  composite  antenna  structure.  The  remainder  of  the 
currents  are  made  up  by  scatter  currents  which  are  produced  by  field  coupling. 
Field  coupling, as  a  first  approximation.  Is  determined  by  the  dominant  current 
distributions,  while  coupling  by  the  scatter  currents  In  general  is  negli¬ 
gible.  Introduction  of  impedances  for  the  characterization  of  structure 
elements  and  their  interaction  permits  utilization  of  network  theory  concepts 
for  the  determination  of  the  junction  currents  and  the  input  impedance  of  the 
antenna.  Formulation  of  all  impedances  by  stationary  expressions  renders  the 
results  insensitive  to  computational  errors  In  the  current  distributions.  As 
demonstrated  by  the  example  given  In  the  paper,  even  rather  crude  approximations 
to  the  dominant  current  distributions  can  yield  good  results. 
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Appendix  1 


Equivalence  between  Current  and  Charge  Excitation 

Consider  a  structure  element  excited  by  an  oscillating  charge  Q  which 
is  placed  at  a  distance  <W3  above  the  (plane)  contact  area  o(Fig.  4).  The 
charge  Q  produces  an  electric  potential  field  -  v$p  which  acts  as  the  pri¬ 
mary  field  for  the  excitation  of  the  structure  element.  The  induced  current 
and  charge  distribution  T,q  radiates  a  Maxwell  field  which  is  characterized 
by  the  retarded  potentials  A  and  $.  The  total  field  satisfies  the  boundary 
condition  *  0: 

[jo»A  +  (v$  +  v$p)]  x  dS  =  0  on  S  and  a 
Current  and  charge  distribution  satisfy  the  continuity  condition 
v*T  +  jwq  *  0  on  S  and  a 
Let 

1  =  1$  on  S,  i  =  on  a,  q  =  q^  on  S,  q  *  qg  on  a  (A1.3) 

A  3  A$  +  Aa  *  =  *s  +  (A1.4) 

where  A^  and  refer  to  the  current  and  charge  distribution  on  S,  and  A^ 
and  +g  to  the  current  and  charge  distribution  on  a.  Since  c«S,  the  contribu¬ 
tion  of  Ag  to  the  total  vector  potential  A  can  be  neglected.  The  charge  distri¬ 
bution  qa  consists  essentially  of  the  counter-charge  to  Q: 

[  q_  da  =  -  Q  (A1.5) 

J  a 


(Al.l) 


(A1.2) 
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There  Is  a  small  additional  induced  charge  on  o  which  is  a  continuation  of 
the  charge  distribution  on  S  into  the  contact  area.  This  charge  can  be  neg¬ 
lected  since  a<<S. 

When  d  approaches  zero,  the  potential  field  of  the  oscillating  charge  Q 
is  compensated  by  that  of  the  counter-charge: 

Thus  ♦  +  io  =  0  ,  7$p  *  7$s  (A1.6) 

This  means,  the  entire  field  is  practically  only  determined  by  the  current  and 
charge  distribution  on  S  which  satisfies  the  boundary  condition 

(jwAs  +  v$s)  x  dS  *  0  on  S  (A1.7) 

and  the  continuity  condition 

v*T^+  jwq g*  0  (A1 .8) 

Moreover,  since  the  net  charge  on  the  structure  element  is  zero,  the  charge  on 
S  is  Q,  and  the  current  flux  through  the  boundary  r  of  the  contact  area  is  ju>Q. 

Thus,  the  current  and  charge  distribution  on  S,  produced  by  the  external 
charge  Q,  are  Identical  with  the  dominant  current  and  charge  distribution  pro¬ 
duced  by  an  Impressed  current  I  *  juQ. 

Since  j«Q  is  the  displacement  current  which  enters  the  structure  element 
at  the  contact  area  it  Is  obvious  that  excitation  by  an  Impressed  displacement 
current  is  equivalent  to  excitation  by  an  impressed  conduction  current. 
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Appendix  2 


Derivation  of  Equation  II_I.j9 

Consider  a  structure  element  with  several  terminals  and  let  k  and  j  be 
any  two  terminals  where  currents  Ik  and  1^  are  impressed.  The  corresponding 
dominant  current  and  charge  distributions  Tk>  qk  and  T^,  q^  produce  the 
fields  Ek  and  Ej  which  satisfy  the  boundary  conditions 


x  dS  =  -  {jwAk  +  7*k)  x  d5  *  0  (A2.1 ) 

Ej  x  dS  =  -  (j<uAj  +  7^)  x  dS  =  5  (A2.2) 

Since  the  currents  ?k  and  are  tangential  to  the  surfaces,  it  follows  from 
the  boundary  conditions 


Tk  dS  *  0 


(A2.3) 


(jwA.  + 


*fc) 


T.  dS  =  0 


(A2.4) 


Using  the  vector  Identity 

v  *  ($?)  *  7$  •  7  +  ♦( v-T)  with  7  •  T  *  -  jwq  (A2.5) 

and  applying  Gauss'  theorem  as  in  (8)  of  Section  II,  equations  (A2.3)  and 
(A2.4)  can  be  written  in  the  form 


C* 


dS  ■  *kk  ‘k 


jw 


L(A> 


T,  + 


♦k  qk)dS 


(A2.6) 
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one  obtains  from  (A2.7)  the  expression  for  Zjk  given  in  the  first  line  of  1 1 1 -1 9. 

The  formulation  in  the  second  line  of  I I I. 19  is  obtained  if  the  potentials 
Ak  and  are  expressed  by  HI.l  and  III. 12  respectively. 
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... 


Let  be  a  dominant  current  distribution  on  the  surface  S1  and 

«sjn  be  the  scatter  current  distribution  on  Sn  produced  by  (n  =  1,  2,  .... 


S*2(r)*£-  l  «Tjn(r')6(r,pI)dSn(rI) 
n=1  J^n 


Multiplying  (A3. 1 )  with  ?n(F)  and  integrating  over  r1 

k 


(A3. 1 ) 


N 

I  *  ^(r)dS1  -£  [  ?J(r)  •  I  [  6T,*n(r*)G(r  ,?')<iSn(?W  Cr) 

Jsi  K  K  w  J$i  K  n*V$n  K 

«^n(r')  •  f  ^(r)G(F,r')dSi(F)dSn(F1) 

n»ljsn  V 

N 

*  I  f  «sTjn(F)  •  A?(r)dSn,  (G(F,F’>  =  G(F\F» 
n=lJsn 

(A3. 2) 


Similarily  it  can  be  shown  that 


I,  »; .  jj  >i «' 


(A3. 3) 


The  proof  for  III. 31  in  the  body  of  the  paper  follows  the  same  outline 


given  above. 


a)  To  prove  that 


Z-  (S*i  +  *<0dS  (A4-1) 

r  Js 

represents  a  stationary  formulation  of  the  intrinsic  impedance,  we  assume  that 
the  dominant  current  distribution  T  has  an  error  of  aT.  The  corresponding  errors 
of  q ,  A  and  $  shal 1  be  denoted  aq ,  aA  and  a$ .  Then 


Z  +  aZ  =  4  (  C(^  +  aA)*(T  +  aT)  +  ($  +  A$)(q  +  aq)]dS  (A4.2) 

r  Js 


The  boundary  condition  for  the  correct  dominant  current  distribution  yields 


(jwA  +  vi)  •  aT  dS  -  0  (A4.3) 

S 


Since  the  dominant  current  distribution  is  the  continuation  of  the  impressed 
current  which  is  assumed  to  be  unchanged,  aT  is  zero  at  the  terminal,  and 
(A4.3)  can  be  written  in  the  form 


jwj  (A  •  aT  +  $aq)dS  *  0 


Using  the  relations 


|  A  •  aT  dS  *  |  aA  ■  T  dS;  |  $aq  dS  *  J  A$q  dS 


(A4.4) 


(A4.5) 
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along  with  (A4.4),  one  obtains 


(A  •  A?  +  $Aq)dS  = 


(aA  •  T  +  A$q)dS 


(A4.6) 


Thus,  from  (A4.1),  (A4.3),  and  (A4.6) 


AZ  ■  f  (aA  *  A1  +  A$Aq)dS 

r  }s 


(A4.7) 


This  means,  AZ  is  of  second  order. 

b)  In  the  case  of  a  mutual  intrinsic  impedance 

zjk '  |s«k  •  !j +  *k‘’j>dS 


(A4.8) 


both  the  dominant  current  distributions  Tk  and  7^  may  have  errors  &7k  and 


Ai..  Thus 

J 


lZjk  *  T^TT  js«lJk  '  *  a*kqj*  *  (Ak  ’  aTj  *  *k  V  *  <*V  lTj  * 

(A4.9) 

Because  the  correct  dominant  current  distributions  satisfy  the  boundary  condi¬ 
tion  E  x  dS  *  5, 

|  (j«Ak  +  7$k)  •  a7 jdS  *  0 


or 


jwj  (Ak  •  a!j  +  ik4qj)dS  *  0 


(A4.10) 
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Furthermore 


Tj  *  1*kt,j)dS  a 


+  |jAqk)dS  =  0 


(A4.ll ) 


From  (A4.10)  and  (A4.ll),  (A4.9)  reduces  to 


Aij  +  A^Aq^dS 


(A4.12) 


Thus  aZj|^  is  of  second  order. 

c)To  prove  that  III. 33  is  a  stationary  expression  for  the  field  coupling 
impedances  we  treat  the  assembly  of  disconnected  structure  elements  like  a 
single  body.  This  means,  when  a  current  is  impressed  on  terminal  (k)  we  con¬ 
sider  the  dominant  current  distribution  ?k  together  with  the  associated  scatter 
currents  which  are  distributed  over  all  the  elements  as  a  dominant  current 
distribution  of  the  system.  The  coupling  impedances  between  any  two  terminals 
can  then  be  formulated  like  mutual  intrinsic  impedances  (A4.8): 


z<f>  ■  *  \  „[<a™  *  *  s<>]ds  (m-'3> 
k,m  IkIm 


The  error  aZ^  produced  by  errors  in  the  current  distributions  T^,  <sTk 


and  ?*,  51*  is  obtained  from  (A4 .12): 


i,* 

AZ(F) 


k,m 


•  Ttif  I  *  iS,m>  +  K  *  ‘♦k’K  *  w,4>lds 

V®  JS" 


(A4.14) 


and  is  of  second  order  This  relation  can  also  be  derived  from  1 1 1. 33 
but  only  in  a  rather  cumbersome  manner. 


Appendix  5 

If  a  current  is  impressed  on  any  terminal  of  a  diakopted  structure  there 
will  be  capacitive  currents  between  the  contact  areas  of  the  disconnected  ele¬ 
ments,  which  have  not  been  considered  in  the  derivation  of  the  field 
coupling  impedances.  One  might  therefore  conclude  that  the  formulas  are  approxi¬ 
mations  which  require  the  gaps  between  adjacent  contact  areas  to  be  so  large 

that  capacitive  currents  are  negligible.  The  purpose  of  this  appendix  is  to 

1 ,i  1  ,i 

show  that  the  expressions  for  Z(F)  and  Z(f)  are  correct  even  if  the  gaps  are 

k,k  k,m 

infinitely  small. 

Figure  25  shows  two  structure  elements,  a  cylindrical  rod  1  and  a  disc  2 

with  the  opposing  contact  areas  a2  and  a-j .  If  a  current  is  impressed  on  the 

terminal  (])  of  the  rod,  there  will  be  a  potential  difference  between  o2  and 
2 

o-|  which.  In  turn,  produces  a  displacement  current  between  these  terminals. 

The  potential  difference  which  is  the  line  integral  of  the  electric  potential 
1  2 

field  between  a2  and  a-j  is  essentially  determined  by  the  charges  on  the  contact 
areas.  If  the  gap  is  made  smaller  and  smaller,  the  potential  difference  ap¬ 
proaches  zero,  and  the  total  current  distribution  becomes  the  dominant  current 
distribution  of  the  interconnected  elements.  As  shown  in  Appendix  1  displace¬ 
ment  currents  at  contact  areas  are  equivalent  to  impressed  currents.  Thus,  the 
situation  discussed  above  is  the  excitation  of  a  diakopted  structure  not  by  one, 
but  by  three  impressed  currents.  To  produce  excitation  by  one  impressed  cur¬ 
rent  in  accordance  with  our  theory  the  displacement  currents  must  be  compen¬ 
sated  so  that  there  is  no  current  flux  from  the  contact  area  onto  the  surface  S 
of  the  element  (S,  by  definition  does  not  c.n  din  the  contact  areas  of  the 
element).  The  magnitude  of  these  compensating  currents  does  not  enter  into 


no 
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Compensation  of  Capacitive  Currents  at  Contact  Areas 


the  analysis  because,  if  the  impressed  currents  of  the  diakopted  structure 
are  identical  with  the  junction  currents  of  the  interconnected  structure  there 
are  no  displacement  currents  between  adjacent  contact  areas  and  the  sum  of 
all  the  compensating  currents  is  zero. 
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